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We investigate the extent to which the uncertainties associated with the propagation of Galactic 
cosmic rays impact upon estimates for the 7-ray flux from the mid-latitude region. We consider 
contributions from both standard astrophysical background (SAB) processes as well as resolved point 
sources. We have found that the uncertainties in the total 7-ray flux from the mid-latitude region 
relating to propagation parameter values consistent with local B/C and 10 Be/ 9 Be data dominate 
by 1-2 orders of magnitude. These uncertainties are reduced to less than an order of magnitude 
when the normalisations of the SAB spectral components are fitted to the corresponding Fermi LAT 
data. We have found that for many propagation parameter configurations (PPCs) our fits improve 
when an extragalactic background (EGB) component is simultaneously fitted to the data. We also 
investigate the improvement in our fits when a flux contribution from neutralino dark matter (DM), 
described by the Minimal Supersymmetric Standard Model, was simultaneously fitted to the data. 
We consider three representative cases of neutralino DM for both Burkert and Einasto DM density 
profiles, in each case simultaneously fitting a boost factor of the DM contribution together with the 
SAB and EGB components. We have found that for several PPCs there are significant improvements 
in our fits, yielding both substantial EGB and DM components, where for a few of these PPCs the 
best-fit EGB component is consistent with recent estimates by the Fermi Collaboration. 
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I. INTRODUCTION 

A plethora of astrophysical data ranging from Galactic 
to cosmological scales indicate the (gravitational) influ- 
ence of otherwise non-interacting particles in our Uni- 
verse. This highly abundant, yet elusive, dark matter 
(DM) , undoubtedly occupies one of the hot seats in cur- 
rent astroparticle physics research. 

This is especially true in light of the recent analyses 
of data from direct detection experiments, which have 
yielded new improved bounds but also intriguing hints of 
a possible signal (see, e.g.,[l], Q). In addition, the Large 
Hadron Collider (LHC) 3], set to probe new areas of 
physics, in particular, supersymmetry (SUSY), which has 
spawned some of the most compelling DM candidates, 
has also been recently activated. (For a review of SUSY 
DM see, e.g., 0H.) 

An alternative strategy to identifying DM is to recon- 
cile discrepancies between current astrophysical predic- 
tions for various radiation (e.g., positrons, antiprotons 
and 7-rays) and current observations, then subsequently 
constrain the properties of DM particles whose annihila- 
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tion (or decay) products may result in such discrepancies. 
To do this, one crucially has to determine how the par- 
ticle and electromagnetic radiation emerging from DM 
annihilations, as well as the radiation emerging from the 
interactions of cosmic rays (CRs) with background nu- 
clei and interstellar radiation fields, propagate through 
the interstellar medium (ISM), whilst specifying the dis- 
tribution of contributing DM and astrophysical sources. 
Unfortunately, the number of Galactic properties rele- 
vant to CR propagation are vast, their effect on these 
indirect signals are largely ambiguous, and the data con- 
straining them (e.g., local nuclear abundance ratios) is 
limited. 

However, the Large Area Telescope (LAT) of the 
recently launched Fermi Gamma-ray Space Telescope 
Q (Fermi LAT henceforth) is designed to observe 7- 
rays with energies 20 McV< E < 300 GeV, which are typ- 
ical of the energies expected from the annihilation prod- 
ucts of the popular group of DM candidates known as 
weakly interacting massive particles (WIMPs) [7] . Fermi 
LAT has acquired data from both the mid-latitude re- 
gion (MLR), 10° < \b\ < 20°^ < I < 360° and 
the Galactic centre (GC) [ll|,[l2[. However, here we only 
utilise data from the MLR since it is likely to be domi- 
nated by 7-rays originating from local sources and hence 
is limited in its contamination from astrophysical sources 
concentrated near the GC 

In this paper we utilise the publicly available GALPR0P 
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package to investigate how different propagation 
parameter configurations (PPCs) affect corresponding 
predictions for (i) the locai abundance of CR nuclei and 
(ii) the 7-ray flux from the MLR. For the former we 
deduce which PPCs are consistent with current measure- 
ments of the abundance ratios B/C and 10 Be/ 9 Be. We 
then utilise these PPCs in the latter analysis, calculating 
the corresponding 7-ray spectra when including a succes- 
sively increasing number of components, including the 
contributions from DM. We assume that the DM solely 
consists of neutralinos described by the Minimal Super- 
symmetric Standard Model (MSSM) (see, e.g., @)), and 
select several benchmark points and density profiles to il- 
lustrate the spread of contributions that can be obtained. 

This paper is organised as follows. In Sec. [II] we briefly 
describe the processes responsible for the main contri- 
butions to the background diffuse 7-ray emission in the 
context of Fermi LAT. In Sec. lIIII we describe the model 
of CR propagation utilised in this study and the rele- 
vant aspects of the GALPROP code used to conduct our 
simulations. In Sec. lIVI we present our results for the 7- 
ray flux emerging from the MLR when considering purely 
astrophysical sources. In Sec.|V] we briefly describe the 
MSSM and the criteria invoked to select our DM candi- 
date points. We then display results for the 7-ray flux 
from the MLR when including contributions from these 
candidate points. Finally in Sec. lVIl we summarise our 
results. 



II. THE BACKGROUND DIFFUSE 
7-RAY EMISSION 

A. Signal components 

The Galactic diffuse 7-ray emission primarily originates 
from the interaction of high energy charged particles 
(known as cosmic rays) with interstellar nuclei and ra- 
diation fields through a variety of physical mechanisms, 
including: bremsstrahlung; inverse Compton scattering 
(ICS); 7r°-decay; and synchrotron emission, each of which 
we now briefly describe. 

Bremsstrahlung occurs in the Galaxy primarily when 
high energy e ± are deflected by the Coulomb field of in- 
terstellar atoms/nuclei, primarily H and He. Here, we 
utilise the original GALPROP code v50.1p to calculate the 
spectra of 7-rays resulting from such interactions. 

Inverse Compton scattering primarily occurs in the 
Galaxy when high energy e interact with photons of the 
interstellar radiation field (ISRF) , up-scattering them to 
higher energies, resulting in energy loss for the travers- 
ing e ± . In this paper, we use a version of GALPROP, that 
we have modified, to calculate the Galactic 7-ray flux 
contribution from the MLR from ICS, taking into ac- 
count three target p hoton distributions constituting part 
of the ISRF [lJ,[II|]: optical (i.e., starlight), the far infra- 
red (FIR) background, and the cosmic microwave back- 



ground (CMB). 

Fluxes of 7-rays from 7r°-decay primarily occur in the 
ISM by pp -chain reactions, generally resulting in the pro- 
duction of two 7-rays. Whilst ICS and bremsstrahlung 
emission are, as we shall see, strongly dependent on the 
propagation parameters describing the diffusion of CRs 
and interstellar nuclei, 7-rays from 7r°-decay, which are 
almost immediately produced where the parent ir° is cre- 
ated, are particularly dependent on the initial distribu- 
tions of nuclei within the Galaxy. Here, we adopt the 
so-called conventional GALPROP model for the initial nu- 
clear abundances 

Finally, synchrotron emission results from the deflec- 
tion of charged particles by the Galactic magnetic field, 
resulting in electro-magnetic radiation emission. Unlike 
bremsstrahlung, ICS and 7r°-decay, the flux contribution 
from synchrotron emission is negligible at the energies 



sensitive to Fermi LAT (see, e.g., [17], [18|)- Hence, we 
ignore such contributions throughout this study. 

In addition to the above, we must account for the con- 
tributions from the hundreds of high energy 7-ray point 
sources identified by Fermi LAT To account for 

these, we utilise results from the analysis presented in 
[l(| (see Sec. |TV] for further details). 

Further, we expect an additional, isotropic, extragalac- 
tic background (EGB) component, which includes the 
sum of contributions from unresolved point sources, dif- 
fuse emission from signatures of large scale structure for- 
mation, diffuse emission resulting from the interaction of 
ultra-high energy CRs and relic photons, etc. (see, e.g., 
[l(J H3, EI| ) . Despite this, if one adopts the estimate of 
the EGB, as recently determined by the Fermi collabo- 
ration [10], one finds that the diffuse 7-ray flux in the 
direction of the GC still dominates the EGB. However, 
the Galactic component decreases rapidly with increas- 
ing latitude to the extent where in the MLR the esti- 
mated EGB becomes comparable at energies E < 1 GeV. 
Therefore, for Fermi LAT 7-ray flux measurements from 
the MLR, observed at energies 20MeV< E <300GeV, if 
one adopts the Fermi estimate the EGB it should be ac- 
knowledged in any corresponding analysis. In Sec. lIVI and 
Sec.|V]we account for the EGB by simultaneously fitting 
a 7-ray flux component, described by the power-law: 
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where A and 7 are fitted parameters, in addition to con- 
tributions from other relevant sources, to the total 7- 
ray flux -E^^f- observed from the MLR by Fermi LAT. 



B. Fermi LAT mid-latitude observations 

The Large Area Telescope on board NASA's Fermi 
Gamma-ray Space Telescope @ has measured the diffuse 
7-ray emission with unprecedented sensitivity and reso- 
lution, being over an order of magnitude more sensitive 
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than its predecessor, the E nerg etic Gamma-Ray Experi- 
ment Telescope (EGRET) |22|. In this study, we utilise 
results from Fermi LAT's observations of the diffuse 7- 
ray emission from the MLR duringFermi LAT's initial 
10 months of scientific observation |8l4ldl|. 

The mid-latitude sky region was selected for the initial 
studies of Fermi LAT since it maximises the contribution 
to the diffuse 7-ray background produced within several 
kiloparsecs of the Sun, hence minimising uncertainties 
associated with CR diffusion. At smaller latitudes, con- 
tributions from sources near to the GC become domi- 
nant, significantly increasing astrophysical ambiguities, 
while at larger Galactic latitudes the emission is increas- 
ingly affected by contamination from charged particles 
misclassified as photons as well as increased uncertain- 
ties in the model used to estimate the background diffuse 
[1, 0| . These results are in strong contention with those 
from EGRET, especially at energies > 1 GeV, where it 
was proclaimed that EGRET observed an excess in dif- 
fuse 7-ray emission when compared to theoretical models 
that correctly reproduce fluxes of directly measured CR 
nucleon and electron spectra [23| . 



III. PROPAGATION OF CR NUCLEONS IN 
THE GALAXY 

A. Description of propagation model 

The composition and energy spectra of the vast majority 
of CRs currently traversing the Galaxy originate from 
the nuclear interactions (also known as spallations) of 
an initial distribution of energetic particles, mainly pro- 
tons and electrons, with nuclei residing in the ISM, and, 
if electrically charged, their electromagnetic interactions 
with the Galactic magnetic field and ISRF. The astro- 
physical sources giving rise to these initial distributions 
of CRs are believed to include mainly supernovae rem- 
nants [24T - I261 ] . and, to a lesser extent, pulsars [27M31I ) . 
compact objects in nearby binary systems [32[ and stel- 
lar winds J3l-[li|. (For a review of potential CR sources 
see, e.g., [36J.) The energetic particles observed by X-ray 
experiments to be accelerating away from these sources 
(see, e.g., [32i|3l| and references therein), propagate sev- 
eral kiloparsecs through the ISM where they lose energy 
via bremsstrahlung and ICS, as discussed in Sec. lII A\ 
as well as ionization and Coulomb interactions, and also 
synchrotron emission as they interact with the Galactic 
magnetic field. These processes modify the energy spec- 
tra and composition of CRs, whilst producing secondary 
particles and 7-ray radiation. 

Throughout the rest of this paper, we refer to the dif- 
fuse 7-ray flux contributions resulting from the interac- 
tion of high energy CRs, originating from the astrophys- 
ical sources mentioned above, with both CR nuclei resid- 
ing in the ISM and the ISRFs as the standard astrophys- 
ical background (SAB) component. 

In order to simulate the propagation of CRs within the 



Galaxy, we utilise the 2D model of the publicly available 
numerical propagation code GALPROP. The intricacies of 
the GALPROP propagation models are described in detail 
elsewhere (see, e.g., [39|, H(|), however in the following 
we briefly summarise the basic features relevant to the 
present study. 

The GALPROP code attempts to numerically solve the 
propagation equation Eq. ([2]) , for a given source distri- 
bution of all CR species, and boundary conditions defin- 
ing the region of propagation, known as the diffusion 
zone, beyond which free-particle escape is assumed. The 
GALPROP propagation equation is written as 
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where ip = ip(r,p,t) is the particle density per unit of 
particle momentum, q(f,p) is the source term, D xx is 
the spatial diffusion coefficient, V is the convection ve- 
locity associated with Galactic winds, re-acceleration is 
described as diffusion in momentum space and is deter- 
mined by the coefficient D ppi Tf is the time-scale for 
nuclear fragmentation, and r r is the time-scale for the 
radioactive decay of nuclei. 

The Galactic magnetic field is generally believed to 
follow the spiral pattern of its stellar population [4l|, [42j , 
and there are several studies that utilise models reflect- 
ing this when investigating CR diffusion (see, e.g., [4l|). 
However, as mentioned in Sec. lII B[ a study of the MLR 
limits the relevant sources to local phenomena and hence 
minimises the effect of the large scale structure of the 
Galaxy. Here we adopt the following model for the Galac- 
tic magnetic field that is default to GALPROP and indepen- 
dent on azimuthal angle: 



B(r,z) = B exp(-r/r - \z\/z ) . 



(3) 



where the scale lengths ro and zq indicate the approxi- 
mate extent of the field. Consequently, it is also appropri- 
ate that the diffusion zone be cylindrically symmetrical, 
and is defined here to be a slab of half-thickness L and 
radius r max . = r = 20kpc [43j]. 

The spatial diffusion coefficient D xx = /3Dq(R/ Ro) a 
is taken here to be a function of particle rigidity R, 
where the relativistic factor j3 — v/c is a consequence of 
a random-walk (diffusion) process, and Rq is a reference 
rigidity set equal to 4 GV. Where diffusive re-acceleration 
is considered, the momentum-space diffusion coefficient 
D pp is taken here to be related to D xx by the following 
relation, derived in 1441. 
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where va is the Alfven velocity, describing the propa- 
gation of low-frequency electromagnetic waves, known as 
Alfven waves, generated by the motion of ions in the ISM 
relative to the Galactic magnetic field [HI, |46[ . 

The convection velocity V(z) incorporated into the 
original GALPROP code to describe the effect of Galac- 
tic winds on CR particle abundances at the periphery 
of the diffusion zone is assumed to increase linearly with 
distance from the Galactic plane (i.e., ^ = ^ = 0, 

> 0), which we adopt here for a range of values of 
. This linear convection model implies a constant adi- 
abatic energy loss, which is consistent with CR-driven 
MHD wind models [H US |1| . 

The GALPROP code incorporates the energy injection 
spectrum for nucleons, which is assumed to be a power- 
law in momentum, dq(p)/dp oc p~" s , with multiple 
breaks: n g = n gi , n g2 below/above a reference rigidity 
with default value R = R. n = 9 GV. This form is shared 
by the injection spectrum for e ± , dq(p) /dp oc p~ e a , again 
with multiple breaks: e g — e go for rigidities R < R\, 
e g = e gi for R\ < R < R2 and e g — e g2 for R > R2, 
where i?i = 4GV and i?i = 10 6 GV are the adopted 
default values. In both cases, the spectra of these parti- 
cles are modified as they lose energy via ionisation and 
Coulomb interactions (nuclei) as well as bremsstrahlung, 
ICS and synchrotron processes (e ± ). 

Lastly, we mention that, when calculating the local 
abundance of CR nuclei, we incorporate the estimated 
effects of the solar cycle on their energies and fluxes fol- 
lowing the prescriptions of the force field approximation 
[49l |. as invoked by Maurin et al. [H(|. Here, the inter- 
stellar energy E IS of a nucleus of charge Ze and atomic 
number A is shifted according to the relation 

E TOA _ E ls _ \Ze\<f> 

where i? TOA is the modulated energy of the nucleus at the 
top of the atmosphere (TO A), and <fi is the solar modu- 
lation potential which varies over the 11 year solar cycle. 
In this study, whilst we assume a time-independent value 
of <f), we utilise a large range of values to illustrate the 
importance of such variations. Given Eq. (jSJ), the rela- 
tionship between the interstellar flux $ IS of interstellar 
nuclei to that at the Earth $ TOA is related by 

frTOA^TOA) / (ffTOA )2 _ (Amp) 2 X 

$/S(£lS) - ^ (£IS)2 _ (Am p y ) ' ^ 

where m p is the proton rest mass. 

B. Evaluation of models: B/C and 10 Be/ 9 Be ratios 

In the previous section we described and highlighted the 
substantial flexibility in the current modelling of the 
propagation of CRs within the Galaxy. In this section 
we illustrate the result of imposing constraints on the val- 
ues of these parameters when invoking the criterion that 



Parameter 


Definition 


Do 


Normalisation of spatial diffusion coefficient 


B 


Galactic magnetic field normalisation 


a 


Spatial diffusion spectral index 


zo 


Scale length of Galactic magnetic field (z-dir.) 


e 30 


Index for e injection spectrum (R < Ri) 


e 31 


Index for e* injection spectrum (Ri < R < R2) 


n si 


Index for nucleon injection spectrum (R < R n ) 


n S2 


index tor nucleon injection spectrum (it > H„) 


VA 


Alfven velocity 


L 


Half-depth of Galactic diffusion zone 


dV 

Ay. 


Convection velocity gradient (z-direction) 


4>o 


Solar modulation potential 



TABLE I: A reference table of the 12 propagation parameters 
involved in our investigation, where rigidities R\ = 4GV, 
R 2 = 10 6 GV and R n = 9 GV. 



PPCs must generate results that are consistent with cur- 
rent measurements of the local interstellar nuclear abun- 
dance ratios B/C and 10 Be/ 9 Be. These ratios are not 
the only two which can be used to potentially constrain 
propagation parameters, but are arguably the most strin- 
gent, which explains their common ap pearance in the 
relevant literature (see, e.g., [171. 139L I5ll458l |). Other pos- 
sible ratios include (Sc+Ti+VJ/Fe, 26 A1/ 27 A1, 36 C1/C1 
and 54 Mn/Mn (see, e.g., [HI, however we neglect 

the use of such data within this study. 

For each of the PPCs that we considered, we utilised 
the GALPROP package to calculate the local abundance ra- 
tios B/C and 10 Be/ 9 Be as a function of nuclear kinetic 
energy, and subsequently compute their corresponding 
Xlar varue (expressed as Xlar P er data point), us- 
ing current measurements of the local abundance ratios 
(LAR) B/C and 10 Be/ 9 Be measurements, defined as 

2 _ g 

where -D^ and o~ij are the respective central values and 
1 a errors of the zth data point of the jth experimental 
data set (either B/C or 10 Be/ 9 Be), IV,- is our predicted 
abundance ratio at the energy corresponding to the data 
point Dij i&ij ■ Hence, the inner sum of the numerator of 
Eq. (0 is over all data points associated with each exper- 
iment, and the outer sum is over all experimental data 
sets. We then normalise this expression by dividing by 
the total number of data points from all experiments con- 
sidered. This quantity is equivalent to a maximum likeli- 
hood method, assuming complete ignorance with respect 
to B/C and 10 Be/ 9 Be measurements, and that the indi- 
vidual measurements are statistically independent from 
one another with Gaussian likelihoods. 

To calculate our results we firstly conducted a broad 
scan of propagation parameter space involving 12 inde- 
pendent parameters. For convenience, we list these pa- 
rameters in TableHl Our broad scan consisted of approx- 
imately 3000 different PPCs, using different values of the 
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Do/10" s 




n 91 


n g2 


5, 6.5, 8.5 


0.2, 0.4, 0.6 


1.7, 2.15, 2.6 


2.0, 2.3, 2.6 


Va 


L 


dV 

dz 




20, 35, 50 


4, 7, 11 


2, 6, 10 


100, 250, 450 



TABLE II: Ranges of the 8 parameters varied in our initial 
broad scan consisting of approximately 3000 PPCs, expressed 
as: Do/10 28 (cm 2 s _1 ), a, n gi , n 92 , va (kms -1 ), L(kpc), 
H (kms" 1 kpc" 1 ) and <f) & (MV). 



8 parameters: D , a, n gi , n g2 , va, L, ^7 and </> Q , dis- 
played in Table|ITJ The values of the 4 remaining pa- 
rameters: B , e ga , e gi and z , were taken to be constants 
equal to the default values specified in the "conventional" 
GALPR0P model (5^G, 1.60, 2.50 and 2kpc, respectively), 
since we found that changing their values from these had 
an insignificant impact on the local nuclear abundance 
ratios. (However, these parameters have a significant ef- 
fect on the 7-ray flux from the MLR, and in Sec. lIVl 
where we present our results for the 7-ray flux, we dis- 
cuss how we utilised PPCs with different values of these 
parameters.) 

Following our initial scan we performed a second higher 
resolution scan, involving a further 3000 points, located 
closely in propagation parameter space to those that pos- 
sessed Xlar < 30 following our initial scan. (We note 
that our high resolution scan on occasion utilised values 
of propagation parameters slightly outside of those ranges 
specified in TablelLTl in order to determine our best- fit re- 
sults. This is illustrated in the Appendix to this paper, 
where we display some of the PPCs utilised in our inves- 
tigation.) 

In Figured] we display the results of our high resolu- 
tion scan. We display the resulting range of B/C (upper 
panel) and 10 Be/ 9 Be ratios (lower panel) in the solar 
neighbourhood, as a function of kinetic energy per nu- 
cleon, and grouped according to their associated values of 
Xlar (innermost to outermost filled regions): Xlar < 3 
(red), 3 < xIar < 4 (blue), 4 < xIar < 10 (green) and 
10 < Xlar < 2 (yellow). In our calcu lati ons of Xlar 
we utilised experimental data from ACE [59( (B/C: black 
circles, 10 Be/ 9 Be: white circles), HEAO-3 [6fl] (B/C: 
white squares), Ulysses [6l| (white diamonds), Voyager 
62] (B/C: white stars, 10 Be/ 9 Be: black stars), the re- 
cent data from the balloon-borne experiment CREAM 
[Si (B/C: black triangles), ISEE-3 [gj ( 10 Be/ 9 Be: white 
triangles), ISOMAX ^ ( 10 Be/ 9 Be: black triangles) and 
other balloon experiments (Balloon) [66l-[68r| ( 10 Be/ 9 Be: 
black squares). Also, in each respective plot we dis- 
play the corresponding spectra associated with the prop- 
agation parameter configuration (PPC) generating the 
smallest value of Xlar resulting from our scan (black 
curves). 

The relatively larger spread in the values of B/C com- 
pared to corresponding results for 10 Be/ 9 Be indicate the 
sensitive nature of the B/C with respect to the selected 
PPC. Despite the purpose of this paper being to high- 
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FIG. 1: B/C (upper panel) and 10 Be/ 9 Be (lower panel) 
in the solar neighbourhood, as a function of kinetic energy 
per nucleon, calculated for various PPCs. Each configura- 
tion is grouped according to its Xlar value per data point 
(innermost to outermost filled regions): Xlar < 3 (red), 
3 < Xlar < 4 (blue), 4 < Xlar < 10 (green) and 10 < 
Xlar < 20 (yellow), calculated using all displayed data sets: 
ACE (B/C: black circles, 10 Be/ 9 Be: white circles), HEAO- 
3 (B/C: white squares), Ulysses (white diamonds), Voyager 
(B/C: white stars, 10 Be/ 9 Be: black stars), CREAM (B/C: 
black triangles), ISEE-3 ( 10 Be/ 9 Be: white triangles), ISO- 
MAX ( 10 Be/ 9 Be: black triangles) and other balloon experi- 
ments (Balloon) ( 10 Be/ 9 Be: black squares). We also display 
the spectra associated with the propagation parameter config- 
uration generating the smallest value of Xlar resulting from 
our scan (thick black curve). 



light the overall uncertainties owing to propagation pa- 
rameters rather than the effects of individual propagation 
parameters, we do mention that, whilst the changes in 
the high energy (i.e., >1 GeV/n) B/C ratio is dominated 
by the effects of the diffusion parameters, the spread in 
the low energy B/C ratio owes to the different values of 
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the solar modulation potential, (j), utilised. This owes 
primarily to the nature of the force field approximation, 
discussed in Sec. lIII| invoked here to estimate the effects 
of the solar wind, where the shift AE = —\Z\tfi/A in en- 
ergy per nucleon, given by Eq. ([5]), is equal to 225 MeV/n 
for both 10 B and 12 C for <f> = 450 MV. Hence, this ex- 
plains the aforementioned importance of <j> at low ener- 
gies, since we expect the effects of solar modulation to 
become significant for energies (per nucleon) near and 
below E ~ \ AE\. 

Compared to the local B/C ratio, the sensitivity of 
10 Be/ 9 Be with respect to PPCs is much less significant, 
indicated in the lower plot of Figure[T] by the much nar- 
rower spread of each Xlar group. Further, we can infer, 
from the closely related spread in 10 Be/ 9 Be of each Xlar 
group at low energies, that this data is much less sig- 
nificant in determining the overall Xlar than the B/C 
data. This clearly reconciles with the poor nature of 
the 10 Be/ 9 Be data relative to the local B/C data. Un- 
fortunately, at higher energies (i.e., > 2GeV), where 
the spread in the 10 Be/ 9 Be spectra from PPCs yielding 
smaller Xlar becomes increasingly narrow, and hence 
would be most significant in constraining propagation pa- 
rameters, there is a complete lack of experimental data. 



IV. 7-RAY FLUX FROM THE MID-LATITUDE 
REGION 

In this section, we present our results for the 7-ray flux 
predicted to emerge from the MLR from astrophysical 
sources and compare them with the corresponding mea- 
surements by Fermi LAT. We utilise a modified version 
of GALPROP to calculate the 7-ray flux for PPCs with 
Xlar < 20, following the high resolution scan of propaga- 
tion parameter space discussed in Sec. lIII Bl We consider 
two different scenarios: 

• Firstly, in Sec. II V Al we consider the 7-ray flux from 
the SAB interactions as well as accounting for con- 
tributions from those point sources identified by 
Fermi LAT in the MLR. 

• Secondly, in Sec. lIVBl we simultaneously fit a 
power-law 7-ray component to the Fermi LAT data, 
intended to represent the EGB flux, in addition to 
the SAB and point source contributions. 



A. SAB and Point Source Contributions 

In Figure^] we display predictions for the total 7-ray flux 
(upper filled regions) from the MLR, consisting of 7- 
ray contributions from: (i) the SAB, (ii) point sources 
identified by Fermi LAT, and (hi) a residual particle 
background consisting of CRs misclassified as 7-rays (re- 
ferred to hereafter as the CR background), as determined 
by the Fermi Collaboration. We compare these pre- 
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FIG. 2: Predicted 7-ray flux from the MLR resulting from 
various astrophysical sources compared with the correspond- 
ing data from Fermi LAT (black data points, with lcr error 
bars). Predictions of the total 7-ray flux (grouped according 
to their Xlar values as in FigureQJ from the combined contri- 
butions from the SAB, point sources identified by Fermi LAT 
(magenta filled region, corresponding to 1 a error range), and 
the Fermi estimated CR background contamination (teal filled 
region). We also display the total 7-ray spectra associated 
with the PPC generating the lowest Xlar value (thick black 
curve), together with its associated spectral components. 



dictions with the corresponding measurements by Fermi 
LAT (black data points, with la error bars). 

To account for point sources, we utilised the Fermi 
LAT point source data presented in [l0|, displayed in 
Figure[2] as the magenta filled region (corresponding to 
the 1 a error range 69] ) . Secondly, we utilised estimates 
for the CR background determined by Monte Carlo sim- 
ulations conducted by the Fermi Collaboration [l(|, dis- 
played in Figure[2] as the teal filled region. 

As mentioned above, to calculate the SAB contribu- 
tion we utilised PPCs with Xlar < 20 following our high 
resolution scan of propagation parameter space. For con- 
venience, the spectra associated with these configurations 
are grouped according to their Xlar va hies using an iden- 
tical colour scheme to that displayed in FigureQ] 

In addition, following the earlier comments made in 
Sec. lIII Bl for each PPC we varied the 4 propagation pa- 
rameters: B ,zo,e go ande 31 , that were found to have no 
discernible influence on the value of Xlar' ano - nence 
were kept constant, but significantly affect the 7-ray flux 
from the MLR. The values of these parameters utilised 
to calculate the 7-ray flux are displayed in Table lllll 

Since the effects of solar modulation are obviously lo- 
cal in nature, leading to the strong ^-dependence of the 
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1.0, 2.0, 3.0 


0.65, 1.60, 1.90 


1.20, 2.50, 3.20 
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TABLE III: Ranges of the 4 propagation parameters, ex- 
pressed as: Bo(jj,G), zq (kpc), e go , e si , that were varied in 
our scan of propagation parameter space when calculating 
the 7-ray flux from the MLR, in addition to those displayed 
in TableHl 



low energy B/C ratio observed in FigureQ] they play no 
role in the spectrum of 7-rays generated by processes 
occurring almost entirely beyond the solar neighbour- 
hood. Therefore, for simplicity, the grouping of the 7- 
ray spectra displayed in Figurc[2] corresponds to the val- 
ues of Xlar associated with a fixed value of <j> = 450 MV. 
To justify this, we investigated how the range of 7- 
ray fluxes changed for the PPCs yielding Xlar < 3 and 
3 < Xlar < 4 , for values 01 = 10 °; 250 and 450 MV. We 
deduced that these changes were minimal in the Fermi 
LAT energy range and have no effect on our conclusions. 

From Figure[2]we can clearly see that the uncertainties 
in the 7-ray flux relating to propagation parameters, i.e., 
approximately 1-2 orders of magnitude in the Fermi en- 
ergy range for our selected PPCs with Xlar < 30, dom- 
inate those associated with other astrophysical sources, 
where in many cases the resulting 7-ray spectra far ex- 
ceeds the Fermi LAT data. This excess can be alleviated 
by adjusting the normalisation of the various spectral 
components of the SAB, which correspond to related ad- 
justments in the normalisation of the e injection spectra 
and initial nuclear abundances. 

For the spectra displayed in Figure[5] the normalisa- 
tions of the spectra and nuclear abundances are all set 
equal to the default values specified in the conventional 
galdef parameter file. This is because such normali- 
sations are irrelevant when fitting the abundance ratios 
B/C and 10 Be/ 9 Be, since the predictions for these ratios 
associated with any given PPC are independent of ab- 
solute fluxes. However, when making predictions for the 
7-ray flux from the MLR, one must take into account 
that the normalisation of the injection spectra of elec- 
trons/positrons as well as the abundances of light nuclei 
can vary. Such variations translate into variations in the 
abundance of supernovae, pulsars, etc. that give rise to 
the initial abundance of electrons/positrons and light el- 
ements within our Galaxy. In this study, we account for 
such uncertainties by allowing for a free variation of the 
three spectral components of the SAB: bremsstrahlung, 
inverse Compton and 7r°-decay emission. As described in 
Sec. Ill Al 7-rays resulting from bremsstrahlung and in- 
verse Compton scattering primarily originate from the 
scattering of electrons/positrons with background nuclei 
or ISRFs. However, 7r°-decay emission primarily origi- 
nates from pions produced in nuclear interactions. 

Consequently, here we represent changes in normalisa- 
tion of the e injection spectra by independently varying 
the combined bremsstrahlung+ICS flux by a factor N, 
relative to the default GALPROP normalisation, and rep- 
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FIG. 3: Best-fit spectra of the total 7-ray flux from the MLR 
for our selected PPCs possessing xlar < 4, following the 
fits described in the text, and grouped according to their 
associated Xbgnd. value: Xbgnd. < 4 (orange), Xbgnd. > 4 
(grey). Each spectrum was fitted to the Fermi LAT data by 
varying the SAB components of each spectra as described in 
the text, using normalisation factors N and M in the range: 
i < N, M < X , for X — 10. The corresponding results when 
using X = 5 were found to be extremely similar (see main 
text). We also display the total 7-ray flux associated with the 
PPC generating the smallest Xbgnd. value (thick black curve), 
together with its corresponding SAB spectral components. 



resent the changes in the initial nuclear abundances by 
varying the normalisation of the 7r°-decay flux by a factor 
M, relative to that corresponding to the results displayed 
in Figure[5J hence, giving rise to the re-normalised SAB 
flux: 
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In FigureEl we display the total 7-ray flux associated 
with PPCs possessing Xlar — 4 a fter renormalising their 
respective spectral components, as described by Eq. 
in order to obtain their best-fits to the Fermi LAT data 
when including the aforementioned contributions from 
point sources and CR background. Best-fit spectra that 
exceeded the ler upper limit of the Fermi LAT data where 
necessarily adjusted to conform with the Fermi LAT data 
whilst maintaining the best-fit possible. The orange and 
grey spectra correspond to PPCs with Xlar — 4 that 
generate best-fit spectra possessing values of the reduced 
X 2 statistic Xb gnd . < 4 and Xb gn d. > 4 respectively, where 
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Xbgnd. is defined as 

s-^N (Dj-Tj) 2 

2 °4 Ml 

Xbgnd. = ~f ■ 9 

n-f 

Here {£).;, cr^} refers to the ith data point of the Fermi 
LAT data, 7$ refers to our corresponding theoretical pre- 
diction, n — 9 is the number of Fermi LAT data points 
utilised, and / is the number of degrees of freedom (see 
below) . 

Here we utilise the reduced % 2 statistic here since in 
what follows we are going to be directly comparing the 
fits of spectra described by different models of varying 
complexity. Consequently, we must invoke a penalty 
for the additional degrees of freedom possessed by each 
model, which we account for here by setting / equal 
to the number of fitting parameters used in addition to 
those utilised to generate the results displayed in Fig- 
ure[2j where no renormalisation took place. Hence, for 
the current scenario, / is to be set equal to the number 
of normalisation factors used (i.e., N and M), and hence, 
is equal to 2. 

We considered two different allowed ranges for the nor- 
malisation factors N and M. In both cases these factors 
were allowed to vary within the range < N, M < X. 
One set of results was generated using X = 10, which cor- 
respond to those displayed in Figure[3J whilst the other 
set of results was generated using X = 5. 

As can be observed from Figure[3l we found that many 
of our selected PPCs provide excellent fits to the Fermi 
LAT data, for both X = 5 and 10, at all relevant ener- 
gies, with Xbgnd va l ues as l° w as 0.653. (For reference, 
TablelVIl in Sec. IVIII lists the values of the best-fit fitting 
parameters, together with their associated x 2 values, for 
the best-fitting PPCs resulting from our described fits.) 

The justification of the exact choice of A is a com- 
plicated question far beyond the scope of this paper, 
depending somewhat on the uncertainties relating to 
the abundance and distribution of various astrophysical 
sources throughout the Galaxy with respect to that mod- 
elled by GALPROP. 

Here we attempt to demonstrate the independence of 
our conclusions between the extreme choice of X = 10 
and the less extreme choice of X = 5, provided that, in 
this scenario, we concern ourselves only with PPCs that 
give rise to a fitted spectra that survive a stringent cut 
of Xbgnd. — 4- Of course, as the value of X increases 
we expect that the number of spectra that survive such 
cuts will increase, since we allow for those spectra with 
a similar shape to the Fermi LAT data, but with un- 
suitable normalisations of its SAB components, to be re- 
normalised to provide the flux necessary to generate a 
small x^gnd value. However, one can observe that be- 
cause all such spectra possess a similar shape (owing pri- 
marily to the power law shape of the Fermi LAT data) the 
spread in the SAB contributions after renormalisation is 
of order 2a (as you would expect for a cut of Xbgnd — 4) 
and hence is not very much different when using either 



X = 10 or X = 5. 

However, as a precaution we will continue to gener- 
ate corresponding results for both X = 10 and X = 5 
throughout, and attempt to explain any significant dif- 
ferences between them and discuss their impact (if any) 
on our conclusions. 



B. SAB + Point Sources + EGB 

As mentioned in Sec. lII A[ if one adopts an extragalac- 
tic background that is similar in magnitude to that es- 
timated by the Fermi Collaboration it is proper to 
incorporate an EGB component in one's attempts to fit 
the Fermi LAT mid-latitude 7-ray data. 

In this section, we present the results of our attempts 
to simultaneously fit the Fermi LAT data with an EGB 
described by the power-law (JTJ 

where, for convenience, E\ ~ 281 MeV is the energy of 
the lowest energy Fermi LAT data point displayed in 
Figure[31 in addition to the SAB and point source com- 
ponents described in Sec. II V Al Once again, we utilised 
those PPCs with Xlar — 4 whose spectra were used to 
perform the fits in Sec. lIV Al (henceforth referred to as 
our "bgnd. only" fits). 

We simultaneously determined the best-fit values of 
the parameters A, 7, N and M for each PPC that gen- 
erated spectra not exceeding the la upper limit of the 
Fermi LAT data, where, once again, we restricted the 
values of N and M to the range < N,M < X, for 
both X = 5 and 10. We constrained our scan to values 
of the EGB normalisation, A, and slope, 7, to the re- 
spective ranges 10 -8 GeVcm~ 2 s _1 sr _1 < A < $1 and 
-4 < 7 < 0, where $i~7x 10 -6 GeVcm~ 2 s _1 sr _1 is 
the la upper limit of the Fermi LAT measurements of the 
total 7-ray flux from the MLR at E 7 — E\. We restricted 
A to values greater than 10~ 8 GeVcm~ 2 s _1 sr _1 since 
we regard smaller best-fit values to indicate a preferen- 
tial absence of an EGB, and restricted A to values smaller 
than $1 since larger values would exceed the Fermi LAT 
data. We restricted 7 to negative values since, irrespec- 
tive of the various astrophysical motivations, it is clear 
from the spectral shape of the Fermi LAT data that only 
negative sloped EGB will benefit the fit to the data, and 
choose a maximum slope of -4 since we believe steeper 
slopes would be difficult to generate from the astrophysi- 
cal sources mentioned in Sec. lII Al that are likely to dom- 
inate the EGB. 

In FigureHJwe display the results of our described fits 
(henceforth referred to as the "Bkg + EGB" fits), for 
X = 10. (We note that there were no significant differ- 
ences between the trends displayed and the correspond- 
ing X = 5 results that would effect our subsequent dis- 
cussion, with the exception that, for X — 5, owing to a 
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FIG. 4: Same as in Figure[3]but now for our bgnd.+EGB fits, 
as described in the text. We also explicitly display the cor- 
responding EGB flux component associated with the best-fit 
spectra possessing the lowest Xdm value (thin dotted curve). 



smaller range for N and M, there are fewer PPCs that do 
not exceed the Fermi LAT data, mostly with Xegb > 4. 
This reaffirms our earlier remarks regarding the indepen- 
dence of our conclusions with respect to the choice of 
X, and hence, for clarity, we henceforth omit displaying 
results corresponding to X — 5.) 

Analogous to the bgnd. only results, the orange and 
grey spectra in Figure|4]correspond to PPCs with Xlar — 
4 that generate best-fit spectra possessing values of the 
reduced x 2 statistic Xegb — 4 and Xegb > 4 respec- 
tively, where Xegb i s identical in form to Xbgnd except 
now Ti includes the EGB flux (JXJ) and / has increased 
to 4 to account for the two additional fitting parame- 
ters A and 7. (For reference, in Tablc lVnl of Scc. lVIII we 
list the values of the best-fit fitting parameters, together 
with their associated x 2 values, for the best-fitting PPCs 
resulting from our described bgnd.+EGB fits.) 

One can clearly observe from Figured that the range 
of fluxes displayed by the orange spectra is much less 
than that displayed in FigureEl This is a consequence 
of the fact that / has increased, increasing the value of 
Xegb f° r those PPCs that fail to improve their fits to 
the data when the additional degrees of freedom associ- 
ated with the EGB are introduced. Hence, those PPCs 
that possess Xegb — 4 still provide reasonable fits to 
the Fermi LAT data when an EGB has been included 
either by: (i) virtue of the EGB (if the EGB component 
is significant) or (ii) because of the excellent fit provided 
by the SAB component irrespective of an EGB. Those 
PPCs that belong to category (i) should be associated 
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FIG. 5: Best-fit values (in addition to several nearby points, 
as described in the text) of the EGB normalisation A and 
slope 7, together with their corresponding value of ^Xegbi f° r 
each of our selected PPCs possessing Xegb < 4 following our 
bgnd.+EGB fits. We also indicate the best-fit EGB estimate 
by the Fermi Collaboration (black star). 



with a significant improvement in their fits to the Fermi 
LAT data when the EGB is included, compared to their 
corresponding bgnd. only fit. Such PPCs should possess 
negative (or very small) values of <5xegb defined by 



r 2 _ 2 2 
"Xegb — Xegb — Xbgnd. 



(10) 



However, those PPCs which reside in category (ii) should 
have positive values of #Xegb owing to the corresponding 
increase of / in Xegb without any significant change in 
the predicted spectrum Tj [7fj| . 

In Figure[5l we display the best-fit values of A and 
7, together with their corresponding value of ^Xegbi f° r 
each of our selected PPCs with Xegb — 4- In order to 
better understand the trends displayed between A, 7 and 
^Xegbj f° r each of the aforementioned best- fit points, in 
Figure[5] we have also plotted several points with values 
of A and 7 close to their corresponding best-fit points, 
again, with their corresponding values of <5xIgb ■ 

We observe that the majority of best-fit values of 
A < 10~ 7 GeVcm~ 2 s~ x sr^ 1 fail to produce an improve- 
ment in the fit to the Fermi LAT data when an additional 
EGB component is included. This owes primarily to the 
size of the error, cr rcs . , associated with the "residual flux" , 
(which we set equal to Di in Eq. © ) equal to the differ- 
ence between the Fermi LAT measurement of the total 
7-ray from the MLR and the sum of the corresponding 
fluxes from point sources and the CR background. Hence, 
we obtain <7 res . by adding the errors of the total, point 
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source and CR background measurements (tJTot. , cps and 
o"CRb respectively) in quadrature: 



cr re 



^PS 



'CRB 

10" (6 " 7) GeVcm 



Tot. 
2 



(11) 



Hence, when one attempts to simultaneously fit the Fermi 
LAT data using the SAB and an associated EGB with 
A <C er rcs . one expects that the improvement in the fit, 
relative to that when just fitting the SAB component 
(i.e., its "bgnd. only" fit), to be marginal. Also, since / 
increases to 4 with our bgnd.+EGB fits we expect that 
PPCs that prefer the EGB with A <C er rcs . to predom- 
inantly have a positive value of <5xegbi which is what 
we observe in Figure[5l Since such EGB components 
have minimal influence in the bgnd.+EGB fits, this also 
explains the lack of any strong correlation between the 
values of A and 7 for A < 10~ 7 GeVcm -2 s _1 sr _1 , in- 
dicating that the quality of the Fermi LAT data is in- 
sufficient to place constraints on components with fluxes 

However, correlations between the best-fit values of A 
and 7 for A > a lcs . allow us to place weak constraints 
on the parameters of the EGB. In Figure[S] we observe 
a preference for slopes with 7 > — 1, which includes 
the best-fit EGB estimate by the Fermi Collaboration: 
i~6x 10" 6 GeVcm^ 2 s _1 sr- 1 and 7 = -0.41 (black 
star) [l(J. Such constraints are not surprising since they 
ultimately result from the slope of the Fermi LAT to- 
tal flux measurements, which one can observe varies ap- 
proximately within this range. Interestingly, the points 
in this region demonstrate a wide range of both positive 
and negative ^XegB' indicating that there are PPCs that 
still provide excellent fits to the Fermi LAT data when a 
significant EGB component is present. 



V. NEUTRALINO DARK MATTER 

In this section we discuss the results obtained from our 
attempts to improve our previous fits to the Fermi LAT 
data by simultaneously fitting a 7-ray flux component 
from SUSY dark matter in addition to those from the 
SAB and EGB. 



A. The Minimal Supersymmetric Standard Model 

Supersymmetry is arguably the current most favoured 
theory of particle physics beyond the Standard Model 
Till . Not only can SUSY solve the hierarchy problem 
72|, in models where R-parity is conserved, the lightest 
neutralino, assumed to be the lightest supersymmetric 
particle (LSP), is a WIMP, whose total annihilation cross 
section often results in a thermal relic density fl h 2 ~ 
0.1, similar to that expected for DM by astrophysical 



observations [73]. Hence.the LSP is, in this respect, an 
excellent DM candidate (3, [1] • 

Because of its large degree of freedom, owing to its 
124 Lagrangian parameters [74j |. and complicated prop- 
erties, it is a common practice to impose within the 
MSSM some well-motivated boundary assumptions, usu- 
ally in the form of grand-unfication conditions on the 
MSSM masses and couplings. The most economical, and 
popular, scenario of this type is the Constrained MSSM 
(CMSSM) [ZjJ. 

On the other hand, even in the general MSSM, only 
some parameters play a real role in determining DM 
properties, including detection rates and fluxes, which 
makes an MSSM-based analysis managable. Further- 
more, the LSP neutralino in the CMSSM (and often in 
other simple unified models) is predominantly gaugino- 
like and, with LEP constraints imposed, heavier than 
roughly 100 GeV, in which case, resulting 7-ray fluxes 
will be uninterestingly small. For these reasons, we con- 
duct our analysis in the framework of the general MSSM. 
On the other hand, it is not our intention here to inves- 
tigate the MSSM in detail, but rather to present typical 
fluxes resulting in the model, along with the uncertain- 
ties associated with the ambiguities in the propagation 
parameters describing the diffusion of CRs. 

Therefore, in this paper we select three representa- 
tive cases of the neutralino in the way described below. 
Firstly, we narrow down the list of MSSM parameters to: 



M 2 , fi, m A 
A T , At, A t 



ttIq, m 



tan/3, 
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(12) 



where M2 denotes the soft wino mass, [i - the 
Higgs/higgsino mass parameter, wla - the pseudoscalar 
Higgs mass, tan/3 - the ratio of the vacuum expectation 
values of the two Higgs fields , mj (I = L, E, Q, U, D) 
- the soft masses of the slepton and squark doublet and 
singlet fields, respectively, and finally Ai (i — r,b,t) - 
the corresponding trilinear couplings. (The remaining 
parameters not specified here are dependent on the val- 
ues of those mentioned above, and related according to 
well-known equations, see, e.g., 72j.) 

In our choice of the above parameters we were guided 
by only varying those which produce the largest effect 
on significant indirect detection parameters, such as the 
thermally-averaged product of the DM annihilation cross 
section and relative velocity, {cT anri v), the yield of 7- 

ray /e^ per annihilation, — d "^ e+ , and the mass of the 
lightest neutralino, Xi\ m x' while simultaneously obey- 
ing the constraints on the DM relic abundance. 

We performed a scan of the above MSSM parameters 
over the ranges as specified in Table lTVl We generated a 
MSSM chain, consisting of 90,000 points, by taking a flat 
prior and peforming a nested-sampling scan using a mod- 
ified version of the publicly available SuperBayeS package 
76]. (An even larger scan over 25 input parameters in 
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< Mi < 2000 


50 < n < 1000 


0<mA< 1000 


2 < tan/3 < 65 


< mj < 2000 


-3000 < Ai < 2000 



TABLE IV: Values, in GeV, of various masses and couplings, 
defined in the main text, used in our scan of MSSM parameter 
space. We take flat prior and likelihood ranges were shown 
on above. 



the general MSSM was performed in [73], although not 
in the context of indirect detection.) We took into con- 
sideration all current collider limits (see Tables 1 and 2 in 
[78j for details) , as well as invoking constraints published 
by the WMAP Collaboration regarding the cosmological 
DM relic abundance I73II. 



B. Point selection 

The 7-ray flux originating from the particle cascade im- 
mediately following each annihilation (which is almost 
entirely due to 7r°-decay), could be considerable when 
observing the MLR. Since 7-rays are not subject to diffu- 
sion processes within the ISM, the differential 7-ray flux 
arriving from a line-of-sight (l.o.s.) at inclination ip rela- 
tive to the direction of the GC is simply given by 
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TABLE V: Table displaying the properties of our three bench- 
mark MSSM DM candidate points. All three points are con- 
sistent with our experimental constraints within the 95% C.L. 
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We see that all the particle physics information re- 
garding DM is conveniently contained within the factor involving for example (a 
(&a.nn.v)m~ 2 -jgr 1 . Hence, for convenience, we define here 
the quantity 



.lGeVHB^^e-r- How- 
ever, with regards to the yields generated, rather 
than invoking a selection criterion for candidate points 
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where 



7,e 



F 1 (E,m x ) is defined as the 



maximum value of E^-^P- within the energy range 
0.1GeV< E 1 < m x , and F e +(E,m x ) is the value of 
-dE e +, where the lower energy limit of 
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0.1 GeV associated with both quantities corresponds to 
the approximate minimum energy of the Fermi LAT sen- 
sitivity range. 

From the above, it is clear that we should expect DM 
candidate points possessing larger values of /susy,7 to 
generate larger 7-ray fluxes within the Fermi sensitiv- 
ity range that originate from the particle cascade fol- 
lowing each annihilation. However, in addition, since a 
large proportion of 7-rays are expected to arise from in- 
teractions involving e ± generated by DM annihilations, 
we expect that the most optimistic 7-ray fluxes will 
be produced by candidate points that generate larger 
yields of e involved in producing 7-rays with energies 
in the Fermi LAT sensitivity range, i.e., candidate points 



.lGcV d_E e+ 

tially selects points with a high yield of and higher 
mass neutralinos, which generally produce a greater num- 
ber of (higher energy) e^. These contribute signifi- 
cantly more to the rate of ICS and bremsstrahlung in 
the Fermi LAT sensitvity range, which, as we shall see, 
dominates the total DM 7-ray flux from the MLR. 

Since we desire our benchmark models to adequately 
reflect the variation of /susy,7 and /suSY,e± 7 as well as 
providing an observable 7-ray flux in the direction of the 
MLR, given the above, we selected the two candidate 
points from our MSSM chain with masses m x ~ 100 GeV 
and 50 GeV corresponding to the points with approxi- 
mately the largest and smallest values of /susy,7 and 
/suSY,e± respectively. For our third benchmark candi- 
date point, we selected the point corresponding to the 
largest posterior probability with respect to all experi- 
mental constraints utilised in generating our MSSM chain 
(denoted hereafter as our best-fit candidate point). 

The properties of our three benchmark candidate 
points are displayed in TablefVl All three points are con- 
sistent with our experimental constraints to within the 
95% C.L. 
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C. Dark matter halo profiles 

To calculate our results we adopt two very different, 
spherically symmetrical halo profiles to describe the DM 
density within the Galactic halo: the Einasto profile (79j 
and the Burkert profile [80]. Owing to indications that 
cuspy profiles are inconsistent with observations, specifi- 
cally regarding the rotation curves of small-scale galaxies 
[8ll - l8q ] , which are more likely to be consistent with den- 
sity profiles possessing flattened cores, we also consider 
the Burkert density profile .80] 
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PBux.( r ) 



[l + (r/r.)] [l + {r/r s ? 



(15) 



Here we use r s = 11.68 kpc and p s = 0.79 GeV cm~ 3 , 
where once again, the latter is determined by normal- 
ising to the local DM density. The values of these pa- 
rameters are consistent with a whole set of dynamical 
informations (see, e.g., [87| and references therein), in- 
cluding, constraints from the trajectories of stars in the 
solar neighbourhood, estimates of the total mass of the 
Galaxy from the motion of satellites in its outer regions, 
and the Galactic rotation curve. The Einasto profile can 
be described by 



pEin.( r ) = Ps ex P 



- 1 



(16) 



where r is the distance from the GC. We use best-fit 
values for the scale radius r s = 21.5 kpc and the slope a = 
0.17, determined from recent numerical simulations of the 
Galactic halo [88[ . Here, we normalise Eq. (fTB) to a local 
density of 0.3GeVcm~ 3 at r = r ~ 8.5 kpc, giving rise 
to a normalisation density factor p s — 0.0538 GeV cm -3 . 



D. Dark matter component 

In this section we present the results of our attempts to 
simultaneously fit a 7-ray flux component from each of 
our DM candidate points, in addition to the SAB and 
EGB, to the Fermi LAT data. 

In Figure[6] we display the residual (i.e. unac- 
counted for) flux following the bgnd.+EGB fits described 
in Sec. lIVBl for PPCs with best-fit spectra possessing 
Xegb — ^ (blue curves). For comparison, we super- 
impose the corresponding un-boosted DM 7-ray fluxes 
(red curves) obtained when using our best-fit (upper), 
gaugino (centre) and mixed (lower) DM candidate points, 
and a Galactic DM density distribution described by our 
Einasto profile (see below for a discussion regarding the 
corresponding results for our Burkert profile). 

We can observe that for energies E 7 < 1 GeV, it is 
inaccurate to approximate the DM 7-ray flux to its ir°- 
decay component, given approximately by Eq. (|13p . which 
is almost entirely independent of propagation parame- 
ters. In each panel of Figure|51 we explicitly illustrate 
the 7-ray flux from 7r°-decays (black dashed curves) . We 
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FIG. 6: Results for the unaccounted for, residual flux (blue 
curves) following the bgnd.+EGB fits described in Sec. lIVBl 
for PPCs with best-fit spectra possessing Xegb __ 4. For 
comparison, we superimpose the corresponding (un-boosted) 
DM 7-ray fluxes (red curves) obtained when using our best- 
fit (upper) , gaugino (centre) and mixed (lower) DM candidate 
points, together with an Einasto DM density profile. In each 
plot we also indicate the 7r°-decay spectral component associ- 
ated with the displayed flux from each respective DM candi- 
date point (dashed black curve). We also display the residual 
flux associated with the PPC whose best-fit spectrum pos- 
sessed the smallest Xegb value following our fits (solid black 
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see that for energies E 1 > 1 GeV, the 7r°-decay flux dom- 
inates the total 7-ray flux, and since this contribution 
is generated local to its source, the associated errors re- 
lating to uncertainties in CR propagation are negligible. 
However, for energies E 1 < lGeV, the 7r°-decay com- 
ponent has virtually no impact on the size or shape of 
the total 7-ray flux from DM, which is now dominated 
by bremsstrahlung and inverse Compton radiation, each 
with uncertainties relating to CR propagation. 

Clearly, one can see that the DM component, whilst 
being consistent with the residual flux at all but a few 
narrow ranges of energies (that can be compensated for 
by appropriately adjusting one or more of our fitting pa- 
rameters) is at least a factor of 10-100 times smaller than 
its corresponding residual flux. Hence, if the 7-ray flux 
from our DM candidate points are to have any significant 
effect when attempting to fit the Fermi LAT data, it is 
clear that such components require substantial enhance- 
ment. To do so, we re-normalised the un-boosted DM 
flux by the commonly adopted boost factor, BF > 1. 

Therefore, for each of the six combinations of our three 
DM candidate points and two DM density profiles, we 
simultaneously determined the best-fit values of the pa- 
rameters N, M, A, 7 and BF, for BF> 1, and using the 
same ranges of values for the other fitting parameters as 
described in Sec. lIVBl 

In Figure [7J we display the results of our fits (hence- 
forth referred to as the "bgnd.+EGB+DM" fits), for 
X = 10, for our mixed DM candidate point when using 
our Burkert (upper panel) , and Einasto (lower panel) DM 
density profiles. (Once again, for the reasons discussed in 
Sec. lIVBl we omit the corresponding results for X = 5.) 
Owing to the flexibility of our fitting procedure and the 
ultimately subdominant effects of the DM component the 
plots corresponding to the other two DM models are simi- 
lar to that displayed, and hence, for clarity, we omit them 
here. Further, we observe that because of these reasons, 
the results displayed, corresponding to the mixed DM 
candidate point for the Burkert and Einasto profiles are 
also extremely similar. 

Analogously with the bgnd. only and bgnd.+EGB fits, 
the orange and grey spectra in FigureUJ correspond to 
PPCs with Xlar — 4 and best-fit spectra possessing val- 
ues of the reduced x 2 statistic Xbm — 4 and Xdm > 4 re- 
spectively, where Xdm * s identical in form to Xegb except 
now Ti includes the DM 7-ray flux and / has increased 
to 5 to account for the boost factor BF. (For reference, 
in TablesmUXIH] of Scc.[VH] we list the values of the 
best-fit fitting parameters, together with their associated 
X 2 values, for the best-fitting PPCs resulting from our 
described bgnd.+EGB+DM fits.) 

Consequently, as described in Sec. lIVBl with regards 
to the bgnd. only and bgnd.+EGB results, because of 
this further increase in / we once again observe a no- 
ticeable reduction in the spread of the best-fit spectra 
with Xdm — 4 compared to that corresponding to the 
bgnd.+EGB results. 

To understand the impact of the DM flux in our 
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FIG. 7: Same as in Figure|4]but now for our bgnd.+EGB+DM 
fits, as described in the text, when using our mixed DM candi- 
date point, together with a Burkert (upper panel) or Einasto 
(lower panel) DM density profile. In each plot we also explic- 
itly display the corresponding DM flux component associated 
with the best-fit spectra possessing the lowest Xdm value (or- 
ange dot-dashed curves). 



bgnd.+EGB+DM fits, and to determine whether the 
presence of a DM component has actually improved the 
fit, one needs to look at the correlation between the pre- 
ferred value of the boost factor BF and the associated im- 
provement in the fit. Analogously to the bgnd.+EGB fits, 
we assess this level of improvement through the quantity 
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FIG. 8: Results for the values of the boost factor BF, EGB 
normalisation A, and corresponding value of Sxdm associated 
with the best-fit points resulting from our bgnd.+EGB+DM 
fits, when using an Einasto DM density profile, together with 
our best-fit (upper left panel), gaugino (upper right panel) 
and mixed (lower panel) DM candidate points. 



S Xdm> defined by 



DM 



Xdm — XeGBj 



(17) 



where negative values of <!>Xdm represent an improvement 
in the reduced x 2 value when the DM contribution is 
included. 

Hence, in Figure we display the best-fit points from 
our bgnd.+EGB+DM fits, with the associated values of 
BF and <!>Xdm indicated, when using our Einasto profile. 
In order to compare the relative importance of the DM 
and EGB components in the best-fit spectra, we also plot 
corresponding values of A on a third (colour scale) axis. 

We can clearly observe that the DM component has a 
less significant effect in improving the bgnd.+EGB fits 
than the EGB did for the bgnd. only fits, with the ma- 
jority of points possessing positive ^XdM' w itli minimum 
<5xdm vahies of -8.95, -6.9 and -8.0 for the mixed, best-fit 
and gaugino models respectively in Einasto Profile. 

We also observe that the distribution of points is very 
similar for each DM model, where in each case there is a 
strong preference for a small range of boost factors ap- 
proximately equal to 40, 1500 and 30 for the best-fit, 
gaugino and mixed DM models respectively. Such corre- 
lations are expected since, analogous to our discussion in 



Sec. lIV B1 regarding the best-fit values of A following our 
bgnd.+EGB fits, these are the approximate boost factors 
necessary to generate respective DM fluxes of order cr rcs ., 
given by Eq. (fTTj) . Hence, we do not expect there to be a 
strong preference for boost factors much less than these 
preferred ranges, with upper limits on BF determined 
by the fact that the total theoretical flux not violate the 
Fermi LAT data. 

The distribution of points corresponding to the use of 
a Burkcrt profile are extremely similar to those displayed 
in Figure[S] Their relationship can be approximately ob- 
tained by rescaling the values of BF displayed in Figure[5] 
by the ratio of the integrals of p 2 x along line of sights, aver- 
aged over the MLR, when using the Burkert and Einasto 
profiles. This ratio is given by 



JMLR d "/l.o., PBur. fcVQttt 
J M LR dn /l.o.s. P%in.( r ^) dl 



0.30, 



(18) 



and is identical to the ratio of the respective DM 7r°-decay 
components when using these two profiles. This relation- 
ship follows because, as can be observed from FigureEJ 
the energies at which the 7r°-decay component dominates 
the DM flux is very similar to those at which the dis- 
played residual flux is calculated. Consequently, since 
the shape of the 7r°-decay spectrum is independent of 
the DM density profile, for a given PPC, the ratio of the 
BF's corresponding to the Einasto and Burkert profiles 
should be approximately equal to the ratio (fT8|) . 

Since the DM annihilation rate is proportional to the 
square of the DM density, a possible source of the above 
enhancements could partially originate from DM sub- 
structures. Recent semi-analytical studies indicate that 
boost factors of up to 20 may arise from substructures 
present in the Galactic halo {89] ■ From Figure[S]we ob- 
serve that the mixed and best-fit candidates together 
with an Einasto profile are the only two of our scenar- 
ios that generate sufficient flux to significantly effect the 
bgnd.+EGB+DM fits when using BF~ 20. One should 
bear in mind however that since it is likely that the distri- 
bution of substructures in the Galactic halo is not congru- 
ent with the smooth DM density profile [88| , the resulting 
boost factor of the DM 7-ray flux will not be spatially 
independent and hence will alter the shape of the DM 
spectrum in a non-trivial way. 

In Figure[51 as we might expect, we observe that points 
with the most negative values of <5xdm have a slight pref- 
erence for larger values of the EGB normalisation where 
the EGB has an increasingly significant effect on the 
bgnd.+EGB+DM fits. In order to further explore such 
correlations, in FigurelH we plot the best-fit EGB param- 
eters corresponding to those points displayed in Figure|8] 

We can clearly observe that the strong correlation 
displayed in Figure[S] between large values of A and 7 
is much diminished, being most distinctive for the re- 
sults corresponding to the gaugino dominated DM point. 
However, interestingly, not only do the best-fit points 
possessing the most negative values of <5xdm possess sig- 
nificant EGB and DM components in order to achieve 
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FIG. 9: Results for the values of the boost factor BF, EGB 
normalisation A, and corresponding value of Sxdm associated 
with the best-fit points displayed in Figure[8] Once again, we 
highlight the location of the Fermi LAT best-fit EGB estimate 
(black stars). 



their fit, their associated EGB component resides in 
the bgnd.+EGB favoured region (which we recall in- 
cludes the Fermi LAT best-fit EGB estimate - indicated 
by the black star in FigureE]), despite the weak cor- 
relation displayed by the bgnd.+EGB+DM fits. (For 
reference, the values {A/10~ 6 , 7, BF, 5xbm\ associated 
with these best-fit points are [1.27,-0.19,55.3,-6.9], 
[1.23,-0.22,1644.5,-8.0] and [1.34,-0.23,48.8,-8.95] 
for the best-fit, gaugino and mixed DM candidate points 
respectively in Einasto profile.) 



VI. SUMMARY 

In this paper we investigated the extent to which the 
uncertainties in the propagation of Galactic CRs impact 
upon our estimates of the 7-ray flux from the MLR. 

For each PPC considered, we calculated the corre- 
sponding Xlar va hie resulting from fits to current exper- 
imental data on the B/C and 10 Be/ 9 Be local abundance 
ratios, and found that a significant number of them were 
consistent with the data, producing \ 2 values as low as 



1.5 per data point. 

We then calculated the corresponding 7-ray spectra 
resulting from standard astrophysical background pro- 
cesses involving energetic CRs and interstellar nuclei or 
ISRFs. We deduced that when the normalisation of the 
various components of the SAB remained fixed, the un- 
certainties in the 7-ray flux generated exceeded several 
orders of magnitude. However, when these normalisa- 
tions were freely adjusted to fit the Fermi LAT data the 
uncertainties in the 7-ray spectra were much reduced, 
with many PPCs producing acceptable fits. 

We then investigated how these fits may be improved 
by simultaneously fitting an EGB component, described 
by a free-form power law, with the SAB to the Fermi 
LAT data. We found substantial improvement in the fits 
of many PPCs, which themselves provided good fits to 
the data, with the largest improvements (i.e. as large 
as <5xegb — —10) being associated with EGB's of order 
A ~ 10~ 7 - 10- 6 GeVcm- 2 s- 1 sr~ 1 at E x = 200 MeV 
and slopes > 7 > — 1, which also includes the EGB 
estimate determined by the Fermi Collaboration. 

We then investigated how these fits may be further im- 
proved by simultaneously fitting a 7-ray flux component 
from dark matter described by three different candidate 
points within the MSSM, together with the SAB and the 
EGB to the Fermi LAT data. We found that the un- 
boosted flux, generated from motivated versions of the 
Burkert and Einasto profiles, was insufficient to impact 
the fits to the data. 

However, when we artificially re-normalised the DM 
flux using boost factors of approximately 40, 1500 and 
30 with our best-fit, gaugino-dominated and mixed 
candidate points together with our Einasto profile (and 
BF's approximately 3.3 times larger when using our 
Burkert profile), we found that there exist many PPCs 
that provide further improvement in their fits to the 
Fermi data. Moreover, there exist PPCs that not only 
provide good fits to the Fermi data when including the 
SAB, EGB and DM components, that are significantly 
better than the corresponding fits when omitting the 
DM flux, but also require both substantial EGB and DM 
fluxes which are crucial in determining the quality of the 
fit with the parameters of the EGB component similar to 
those estimated previously by the Fermi Collaboration. 
A possible source of these enhancements could arise 
from local DM substructures within the Galactic halo, 
thought to be able to give rise to boost factors as large 
as 20. 
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TABLE VI: Tabulated values of those PPCs utilised in our bgnd. only fits (i.e. Xlar) that generated the smallest values of 
Xbgnd. f° r X = 10. We also provide the best-fit values of the normalisation parameters N and M of each PPC as well as the 
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0.4 


1 


0.65 


3.6 


1.7 


2 


35 


11 


6 


450 


3.94 


0.603 


2.2 


-8.92 


-3.39 


1.4 


0.399 


5 


5 


0.38 


2 


1.6 


2.5 


1.79 


2.25 


35 


4 


2.5 


450 


1.64 


1.02 


1.25 


-5.75 


-0.354 


1.42 


-1.57 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2 


20 


7 


2 


100 


2.82 


2.37 


1.34 


-9.82 


-3.65 


1.42 


0.407 


5.2 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2 


450 


1.72 


0.792 


1.24 


-5.72 


-0.343 


1.44 


-2.12 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.15 


35 


4 


2.5 


450 


1.63 


1.02 


2.25 


-6.06 


-1.65 


1.44 


-0.234 


6.5 


2.5 


0.4 


3 


0.65 


3.6 


1.7 


2.3 


35 


11 


10 


450 


3.35 


0.641 


1.36 


-6.39 


-0.0441 


1.44 


-1.8 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


6 


250 


3.46 


1.44 


1.8 


-5.66 


-0.713 


1.45 


-2.64 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2.3 


35 


4 


6 


450 


3.76 


0.562 


1.17 


-6.04 


-0.163 


1.46 


-5.76 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.27 


35 


4 


2.5 


450 


1.63 


0.599 


1.28 


-5.62 


-0.387 


1.46 


-1.96 


5 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


2.6 


0.95 


1.02 


-5.73 


-0.326 


1.46 


-2.82 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


42 


14 


8 


450 


2.79 


0.903 


1.07 


-5.91 


-0.255 


1.49 


-2.67 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


4 


450 


1.9 


0.868 


0.914 


-5.69 


-0.337 


1.5 


-2.3 



TABLE VII: Same as for Table lVTl but now for our bgnd.+EGB fits. For each PPC listed we also provide the best-fit values of 
the EGB normalisation A (GeVcm -2 s _1 sr" 1 ) and slope 7, as well as the corresponding values of Xegb an d <5xegb- 
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Propagation Parameter Combination 



Burkert Profile 



D 


Bo 


a 


80 


e 90 


e si 


n 31 


n 32 


VA 


L 


dV 

dz 


4>@ 


i 

Xlar 


N 


M 


logio(A) 


7 


logio(BF) 


'2 

Xdm 


Sxdm 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


4 


2 


450 


1.98 


0.672 


1.94 


-7.24 


-0.0566 


1.94 


0.842 


-1.36 


5 


2.5 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.28 


1.3 


-9.24 


-2.59 


2.15 


0.979 


-2.03 


5.75 


10 


0.34 


1 


0.65 


3.6 


1.82 


2.36 


36 


4 


2.5 


450 


3.17 


0.552 


0.663 


-5.96 


-0.188 


2.12 


1.02 


-1.37 


5.75 


10 


0.34 


3 


0.65 


3.6 


1.82 


2.36 


38 


4 


3 


450 


2.92 


0.556 


0.596 


-5.99 


-0.171 


2.18 


1.07 


-1.87 


5 


2.5 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.33 


1.12 


-7.92 


-1.04 


2.22 


1.1 


-1.9 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


6 


250 


3.8 


1.42 


2.04 


-6.07 


-0.954 


1.07 


1.11 


0.156 


5.75 


2.5 


0.35 


3 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


3.77 


0.571 


0.814 


-5.86 


-0.256 


2.01 


1.11 


-1.19 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


6 


250 


1.71 


1.73 


1.91 


-6.61 


-1.83 


1.06 


1.14 


0.243 


8 


5 


0.4 


2 


1.6 


2.5 


1.8 


2.36 


40 


14 


8 


450 


3.24 


0.84 


1.22 


-5.69 


-0.438 


1.76 


1.16 


0.0817 


5 


10 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.04 


1.12 


-8.05 


-3.55 


2.36 


1.17 


-3.69 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


11 


10 


450 


3.35 


2.13 


1.51 


-6.22 


-1.73 


0.884 


1.18 


0.228 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


40 


14 


8 


450 


3.34 


1.4 


1.88 


-6.4 


-2.46 


1.29 


1.19 


0.095 


6.5 


10 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


11 


10 


450 


3.31 


0.648 


1.64 


-7.65 


-0.869 


2.21 


1.23 


-2.43 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


2 


250 


3.03 


2.25 


1.36 


-7.29 


-2.02 


1.3 


1.28 


0.107 


5 


5 


0.36 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2.5 


450 


1.8 


0.9 


1.23 


-5.74 


-0.359 


1.59 


1.34 


-0.303 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


2 


250 


2.15 


1.31 


1.72 


-6.21 


-0.426 


1.39 


1.35 


0.398 


5 


2.5 


0.4 


3 


0.65 


3.6 


2.15 


2 


35 


4 


10 


450 


3.8 


0.735 


1.53 


-6.98 


-0.994 


2.25 


1.38 


-1.14 


4.8 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


4 


2 


450 


1.87 


0.708 


1.08 


-5.62 


-0.396 


1.75 


1.38 


-0.587 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


20 


7 


2 


250 


3.28 


1.71 


1.82 


-5.9 


-1.21 


0.892 


1.4 


0.575 


5 


2.5 


0.4 


1 


0.65 


3.6 


2.15 


2 


35 


4 


10 


450 


3.8 


0.752 


1.55 


-7.55 


-1.27 


2.21 


1.42 


-1.75 


7 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


40 


14 


8 


450 


2.41 


0.74 


0.933 


-5.72 


-0.349 


1.85 


1.43 


-0.402 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


34 


14 


3 


450 


2.59 


0.572 


1.09 


-5.67 


-0.385 


1.74 


1.43 


-0.755 


5 


10 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.05 


1.04 


-9.41 


-3.54 


2.35 


1.43 


-2.58 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


6 


250 


2.79 


2.25 


1.12 


-8.57 


-1.83 


1.68 


1.45 


0.549 


6.5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


2 


100 


3.91 


0.719 


1.15 


-5.7 


-0.356 


1.16 


1.45 


0.065 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.6 


20 


7 


6 


250 


1.78 


0.545 


0.935 


-5.76 


-0.322 


1.99 


1.45 


-0.509 


5.75 


5 


0.35 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


3.77 


0.904 


1.31 


-5.79 


-0.38 


1.74 


1.46 


-0.614 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


42 


14 


8 


450 


2.79 


0.952 


0.858 


-5.76 


-0.341 


1.8 


1.47 


-0.278 


5 


2.5 


0.4 


3 


0.65 


3.6 


1.7 


2 


35 


4 


6 


450 


3.97 


0.916 


1.3 


-6.59 


-0.0576 


1.66 


1.47 


-0.106 


8 


10 


0.45 


1 


0.65 


3.6 


1.8 


2.36 


40 


14 


5.5 


450 


2.13 


0.386 


0.64 


-6.22 


-0.13 


2.33 


1.48 


-5.67 


5.1 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


2.5 


450 


1.66 


1.04 


1.36 


-5.84 


-0.344 


1.48 


1.49 


-0.128 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


35 


14 


3 


450 


2.68 


1.25 


1.62 


-6.44 


-0.232 


1.26 


1.5 


0.371 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.15 


35 


4 


2.5 


450 


1.63 


1.24 


2.12 


-6.39 


-3.48 


1.45 


1.5 


0.0606 


5.2 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2.5 


450 


1.83 


0.798 


1.19 


-5.68 


-0.386 


1.59 


1.54 


-0.299 


6.5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2 


35 


11 


6 


250 


2.7 


0.909 


2.19 


-6.28 


-2.39 


1.5 


1.55 


-0.011 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


7 


10 


450 


3.78 


1.69 


1.45 


-6.36 


-2.22 


1.9 


1.58 


-0.412 


8.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


35 


11 


2 


450 


2.58 


1.07 


0.64 


-6.07 


-0.215 


2.06 


1.59 


-0.716 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


40 


4 


2 


450 


2.01 


0.653 


1.05 


-5.79 


-0.341 


2.01 


1.6 


-0.608 


4.8 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2.5 


450 


2.19 


0.63 


1.09 


-5.62 


-0.404 


1.8 


1.61 


-0.478 


7.5 


10 


0.45 


3 


0.65 


3.6 


1.8 


2.36 


40 


14 


8 


450 


2.18 


0.451 


0.398 


-5.97 


-0.162 


2.28 


1.63 


-1.61 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.35 


35 


4 


2.5 


450 


1.65 


0.705 


0.981 


-5.64 


-0.391 


1.88 


1.63 


-0.183 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


7 


6 


450 


1.92 


1.83 


1.33 


-6.55 


-3.39 


1.86 


1.65 


0.704 


5 


5 


0.4 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


3 


450 


1.82 


0.712 


1.21 


-5.76 


-0.385 


1.94 


1.66 


-1.35 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.18 


36 


4 


3 


450 


3.69 


3.15 


1 


-7.98 


-1.42 


0.598 


1.67 


0.416 


8 


10 


0.45 


3 


0.65 


3.6 


1.8 


2.36 


35 


10 


3 


450 


1.76 


0.402 


0.618 


-6.2 


-0.109 


2.31 


1.67 


-5.25 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


7 


6 


250 


2.27 


2.76 


0.798 


-8.7 


-0.862 


1.87 


1.68 


-0.747 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2 


35 


11 


10 


250 


1.62 


1.2 


2.06 


-6.63 


-2.98 


1.31 


1.68 


0.248 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.1 


35 


4 


2 


450 


1.92 


1.95 


1.66 


-6.57 


-3.13 


1.29 


1.68 


0.228 


8.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


11 


2 


450 


2.08 


1.34 


1.55 


-6.39 


-3.41 


1.89 


1.69 


0.17 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


34.5 


4 


2.5 


450 


1.67 


0.816 


1.52 


-5.8 


-0.471 


1.84 


1.72 


-0.258 



TABLE VIII: Same as for Table lVIII but now for our bgnd.+EGB+DM fits when using our best-fit DM candidate point with 
our Burkert DM density profile. For each PPC listed we also provide the best-fit values of the best-fit boost factors BF as well 
as the corresponding values of Xdm an< i ^Xdm- 
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Propagation Parameter Combination 



Einasto Profile 



D 


Bo 


a 


«o 


e 90 


e 91 


n 91 


n 92 


VA 


L 


dV 

Hz 


<t>Q 


'2 

Xlar 


N 


M 


logio(^) 


7 


logio(BF) 


'2 

Xdm 


<5Xe>m 


6.5 


2.5 


0.4 


3 


0.65 


3.6 


1.7 


2 


35 


7 


6 


450 


1.93 


0.71 


1.85 


-10 


-1.07 


1.52 


0.492 


-2.08 


5.75 


2.5 


0.34 


3 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


3.52 


0.598 


0.893 


-5.97 


-0.226 


1.55 


0.86 


-0.749 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


2 


250 


2.62 


1.37 


0.842 


-5.75 


-0.361 


0.931 


0.953 


-0.246 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


6 


250 


3.64 


2.11 


0.899 


-6.02 


-0.351 


0.998 


0.965 


-0.0893 


5 


2.5 


0.34 


3 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


2.6 


0.545 


0.804 


-6.03 


-0.211 


1.66 


0.989 


-1.5 


5 


2.5 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.27 


1.3 


-7.64 


-0.853 


1.6 


1.03 


-1.98 


8 


2.5 


0.45 


3 


0.65 


3.6 


1.8 


2.36 


35 


14 


2.5 


450 


3.81 


0.403 


0.687 


-6.07 


-0.144 


1.65 


1.04 


-2.58 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


3.52 


0.899 


1.29 


-5.71 


-0.379 


0.983 


1.13 


-0.232 


6.5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


2 


100 


3.91 


1.38 


1.14 


-6.13 


-0.299 


1.14 


1.14 


-0.246 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


2 


250 


2.15 


1.2 


1.76 


-6.2 


-0.302 


0.47 


1.14 


0.195 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


6 


250 


3.8 


1.06 


2.28 


-6 


-1.12 


0.663 


1.14 


0.191 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


6 


250 


2.82 


0.848 


0.738 


-5.84 


-0.298 


1.46 


1.16 


-1.01 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


40 


14 


8.5 


450 


3.97 


0.785 


1.15 


-5.73 


-0.354 


0.97 


1.16 


-0.252 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


6 


250 


1.71 


1.61 


2 


-6.66 


-3.5 


0.225 


1.17 


0.276 


5 


10 


0.4 


1 


0.65 


3.6 


1.7 


2.3 


35 


7 


10 


450 


3.78 


0.569 


0.583 


-6.21 


-0.0958 


1.7 


1.18 


-1.52 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


5 


450 


1.83 


0.601 


1.09 


-5.62 


-0.367 


0.799 


1.19 


-0.0714 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


40 


14 


8 


450 


3.34 


1.2 


1.96 


-6.09 


-1.38 


0.931 


1.19 


0.0976 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


2.5 


450 


3.17 


0.828 


1.21 


-5.65 


-0.377 


0.773 


1.2 


-0.0317 


5 


2.5 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.35 


1.03 


-8.85 


-3.65 


1.7 


1.23 


-1.78 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


35 


14 


3 


450 


2.68 


0.929 


1.5 


-5.8 


-0.421 


0.652 


1.25 


0.122 


6.5 


10 


0.4 


3 


0.65 


3.6 


2.15 


2 


35 


11 


10 


450 


1.9 


0.528 


1.77 


-7.33 


-0.0875 


1.55 


1.26 


-1.36 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


2 


250 


3.03 


1.74 


1.69 


-6.82 


-3.7 


0.871 


1.29 


0.119 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


11 


10 


450 


3.35 


2.7 


1.1 


-6.86 


-2.95 


0.276 


1.33 


0.443 


5 


2.5 


0.4 


1 


0.65 


3.6 


2.15 


2 


35 


4 


10 


450 


3.8 


0.764 


1.5 


-7.95 


-1.33 


1.69 


1.36 


-1.81 


5.2 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2.5 


450 


1.83 


0.668 


1.37 


-5.74 


-0.358 


0.963 


1.38 


-0.185 


5 


2.5 


0.4 


3 


0.65 


3.6 


2.15 


2 


35 


4 


10 


450 


3.8 


0.676 


1.73 


-7.12 


-2.15 


1.58 


1.38 


-1.14 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


2.4 


450 


1.65 


0.783 


1.4 


-5.74 


-0.39 


1.1 


1.39 


-0.418 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.27 


35 


4 


2.5 


450 


1.63 


0.599 


1.28 


-5.62 


-0.387 


0.835 


1.39 


-0.066 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


7 


10 


450 


3.78 


1.06 


1.46 


-5.77 


-0.516 


1.15 


1.4 


-0.59 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2.3 


35 


4 


2 


450 


1.86 


0.697 


0.545 


-6.17 


-0.088 


1.55 


1.41 


-0.899 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


34 


4 


2 


450 


1.71 


0.816 


1.28 


-5.78 


-0.318 


0.867 


1.45 


-0.0798 


5 


10 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.05 


1.06 


-7.71 


-3.44 


1.82 


1.46 


-3.39 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


3.5 


450 


1.69 


0.872 


1.57 


-5.88 


-0.376 


1.1 


1.47 


-0.4 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


6 


250 


2.79 


2.16 


1.17 


-7.62 


-0.688 


1.14 


1.48 


0.579 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


5.5 


450 


2.14 


1.57 


0.934 


-6.15 


-0.27 


1.2 


1.48 


0.391 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


3 


450 


1.65 


0.992 


1.33 


-5.78 


-0.383 


1.03 


1.48 


-0.302 


8 


5 


0.45 


2 


1.6 


2.5 


1.7 


2.36 


35 


14 


3 


450 


2.41 


1.26 


0.992 


-6.13 


-0.247 


1.27 


1.49 


-0.333 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


7 


10 


450 


3.21 


0.683 


2.04 


-6.69 


-2.4 


0.841 


1.5 


-0.228 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


40 


14 


5.5 


450 


2.13 


0.619 


0.993 


-5.75 


-0.371 


1.43 


1.53 


-0.867 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


4 


2 


450 


1.86 


0.766 


1.19 


-5.69 


-0.407 


1.28 


1.53 


-0.175 


8 


2.5 


0.45 


3 


0.65 


3.6 


1.8 


2.36 


40 


14 


5.5 


450 


2.13 


0.398 


0.571 


-5.9 


-0.225 


1.6 


1.53 


-1.48 


5 


10 


0.4 


1 


0.65 


3.6 


1.7 


2.3 


35 


4 


6 


450 


3.76 


0.598 


0.748 


-6.13 


-0.215 


1.76 


1.56 


-0.766 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.15 


35 


4 


2.5 


450 


1.63 


1.08 


2.02 


-5.95 


-1.13 


1.12 


1.58 


0.135 


8 


2.5 


0.45 


1 


0.65 


3.6 


1.8 


2.36 


35 


14 


3 


450 


2.93 


0.404 


0.597 


-5.96 


-0.183 


1.62 


1.59 


-1.3 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.35 


35 


4 


2 


450 


1.71 


0.715 


0.99 


-5.68 


-0.344 


1.2 


1.59 


-0.349 


5 


10 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.04 


0.715 


-6.42 


-0.167 


1.73 


1.61 


-2.41 


6.5 


2.5 


0.4 


1 


1.9 


2.95 


1.7 


2.6 


35 


11 


10 


450 


3.49 


0.54 


0.732 


-5.84 


-0.267 


1.56 


1.61 


-2.64 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


35 


7 


6 


450 


2.02 


0.985 


0.739 


-6.01 


-0.231 


1.47 


1.61 


-0.962 


5 


10 


0.6 


1 


0.65 


3.6 


1.7 


2.6 


35 


11 


10 


250 


1.5 


0.293 


0.325 


-5.97 


-0.169 


1.75 


1.62 


-5.71 


6.5 


2.5 


0.4 


3 


0.65 


3.6 


1.7 


2 


35 


11 


10 


450 


3.31 


0.793 


1.86 


-7.84 


-3.8 


1.08 


1.63 


-0.673 



TABLE IX: Same as for Table lVIHl but now using our Einasto DM density profile. 



22 



Propagation Parameter Combination 



Burkert Profile 



Do 


B 


a 


Zo 


6 90 


e 91 


n 31 


n 92 


VA 


L 


dV 

d-/. 


4>B 


2 

Xlah. 


N 


M 


logio(A) 


7 


logio(BF) 


'2 

Xdm 


<5Xdm 


5 


2.5 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.24 


1.3 


-7.88 


-0.826 


3.61 


0.743 


-2.27 


5.5 


10 


0.34 


3 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


2.25 


0.532 


0.621 


-6.03 


-0.169 


3.71 


0.795 


-1.62 


6.5 


10 


0.4 


3 


0.65 


3.6 


1.7 


2.6 


35 


11 


10 


450 


3.49 


0.497 


0.414 


-6 


-0.165 


3.72 


0.917 


-2.64 


5 


2.5 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.26 


1.31 


-7.33 


-1.08 


3.52 


0.957 


-2.05 


8 


2.5 


0.45 


1 


0.65 


3.6 


1.8 


2.36 


40 


14 


5.5 


450 


2.13 


0.403 


0.53 


-5.95 


-0.166 


3.6 


0.969 


-2.32 


5 


10 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1 


1.27 


-9.34 


-3.03 


3.65 


0.973 


-3.04 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


40 


14 


8.5 


450 


3.97 


0.697 


1.12 


-5.67 


-0.365 


2.94 


1.02 


-0.322 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


4 


6 


450 


3.97 


0.805 


1.95 


-8.23 


-2.17 


3.23 


1.04 


-0.347 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


6 


250 


3.8 


1.17 


2.26 


-6.07 


-1.32 


2.55 


1.07 


0.119 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


7 


6 


450 


1.92 


0.949 


1.38 


-5.7 


-0.455 


2.98 


1.08 


0.142 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


6 


250 


1.71 


2.2 


1.66 


-7.33 


-2.71 


1.32 


1.1 


0.21 


5 


2.5 


0.4 


3 


0.65 


3.6 


1.7 


2 


35 


4 


6 


450 


3.97 


0.847 


1.74 


-6.85 


-0.0873 


2.96 


1.11 


-0.473 


5 


10 


0.4 


1 


0.65 


3.6 


1.7 


2.3 


35 


7 


10 


450 


3.78 


0.56 


0.313 


-6.06 


-0.127 


3.71 


1.12 


-1.58 


8.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


35 


11 


2 


450 


2.58 


0.489 


0.627 


-5.57 


-0.352 


3.2 


1.14 


-1.06 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


2 


250 


2.15 


1.22 


1.68 


-6.11 


-0.296 


1.8 


1.14 


0.194 


5.75 


10 


0.34 


1 


0.65 


3.6 


1.82 


2.36 


36 


4 


2.5 


450 


3.17 


0.533 


0.444 


-5.85 


-0.241 


3.7 


1.17 


-1.22 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


11 


10 


450 


3.35 


2.79 


1.07 


-7.08 


-2.8 


1.76 


1.17 


0.287 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


40 


14 


8 


450 


3.34 


1.23 


2.03 


-6.3 


-3.07 


2.64 


1.18 


0.0997 


5 


2.5 


0.34 


1 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


2.6 


0.584 


0.504 


-5.85 


-0.264 


3.68 


1.19 


-1.68 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


7 


2 


250 


3.2 


0.646 


0.97 


-5.68 


-0.36 


3.3 


1.2 


0.0573 


8 


2.5 


0.45 


3 


0.65 


3.6 


2 


2.36 


40 


14 


8 


450 


2.07 


0.395 


0.647 


-6.02 


-0.137 


3.53 


1.22 


-2.67 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


2 


250 


3.03 


1.91 


1.55 


-6.86 


-1.94 


2.95 


1.27 


0.0904 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


5.5 


450 


2.14 


1.48 


1.02 


-6.31 


-0.187 


3.23 


1.28 


-1.09 


8.5 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


3 


450 


1.8 


0.653 


0.955 


-5.68 


-0.372 


3.35 


1.29 


-0.911 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


3 


450 


2.93 


0.579 


0.876 


-5.6 


-0.373 


3.22 


1.3 


-0.374 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


20 


7 


6 


250 


2.21 


1.62 


1.83 


-5.73 


-0.793 


1.81 


1.3 


0.246 


5 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


2.6 


0.804 


1.01 


-5.68 


-0.349 


3.12 


1.31 


-0.148 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


4 


2 


450 


1.98 


0.807 


0.961 


-6.31 


-0.151 


3.48 


1.32 


-0.876 


8 


10 


0.45 


3 


0.65 


3.6 


1.8 


2.36 


35 


14 


5 


450 


1.83 


0.441 


0.511 


-6.03 


-0.158 


3.67 


1.32 


-2.58 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


10 


3 


450 


1.76 


0.609 


1.04 


-5.78 


-0.312 


3.25 


1.32 


-1.03 


5 


5 


0.6 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


6 


250 


2.04 


0.923 


1.13 


-5.97 


-0.261 


3.27 


1.34 


-0.599 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.1 


35 


4 


2 


450 


1.92 


1.24 


1.98 


-5.86 


-0.794 


2.75 


1.37 


-0.0904 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


20 


7 


2 


250 


3.28 


1.01 


2.37 


-5.83 


-1.42 


2.24 


1.37 


0.543 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


2.5 


450 


3.17 


1.3 


1.03 


-5.75 


-0.412 


3.27 


1.37 


-0.0017 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


3.5 


450 


1.69 


0.938 


1.46 


-5.84 


-0.398 


3.16 


1.41 


-0.814 


8 


10 


0.45 


1 


0.65 


3.6 


1.8 


2.36 


35 


14 


5.5 


450 


2.14 


0.481 


0.326 


-6 


-0.13 


3.65 


1.45 


-1.89 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


3 


450 


1.65 


0.896 


1.37 


-5.77 


-0.391 


3.03 


1.46 


-0.325 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


3.52 


1.06 


1.2 


-5.76 


-0.449 


3.37 


1.46 


0.0961 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


6 


250 


2.79 


2.44 


0.944 


-9.65 


-2.01 


3.19 


1.46 


0.231 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.15 


35 


4 


2.5 


450 


1.63 


1.09 


2.08 


-6.12 


-1.33 


3.12 


1.48 


-0.0124 


8 


5 


0.5 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


3 


250 


2.32 


0.559 


0.885 


-5.77 


-0.327 


3.47 


1.48 


-0.904 


6.5 


10 


0.4 


1 


0.65 


3.6 


1.7 


2.3 


35 


11 


10 


450 


3.35 


0.507 


0.828 


-5.96 


-0.207 


3.41 


1.49 


-1.27 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


4 


2 


450 


1.86 


0.91 


0.982 


-5.65 


-0.354 


2.91 


1.5 


-0.288 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


20 


7 


2 


100 


2.75 


1.14 


1.13 


-6.36 


-0.173 


3.46 


1.5 


-1.49 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


2 


250 


2.62 


2.21 


0.691 


-6.22 


-0.386 


3.33 


1.53 


0.331 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


35 


14 


3 


450 


2.68 


1.01 


1.88 


-6.2 


-2.48 


3.26 


1.53 


0.408 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


6 


250 


2.82 


0.911 


0.775 


-6.09 


-0.198 


3.48 


1.53 


-0.628 


8 


2.5 


0.45 


1 


0.65 


3.6 


1.8 


2.54 


35 


14 


3 


450 


3.23 


0.326 


0.437 


-5.82 


-0.248 


3.6 


1.54 


-3.45 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2 


35 


11 


10 


250 


1.62 


0.994 


2.07 


-6.09 


-0.72 


2.8 


1.54 


0.108 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


2.4 


450 


1.65 


0.854 


1.34 


-5.76 


-0.39 


3.01 


1.56 


-0.243 



TABLE X: Same as 



for Table lVlTIl but now for 



our gaugino-' 



dominated DM candidate point. 



2.3 



Propagation Parameter Combination 



Einasto Profile 



Do 


B 


a 


Zo 


e 90 


e 91 


n 31 


n 92 


VA 


L 


dV 

d-/. 


4>B 


■2 

Xlar 


N 


M 


logio(A) 


7 


logio(BF) 


'2 

Xdm 


<5Xdm 


6.5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


11 


6 


450 


3.94 


0.603 


2.2 


-8.92 


-3.39 


2.68 


0.45 


-0.949 


5 


2.5 


0.34 


1 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


2.6 


0.608 


0.565 


-5.97 


-0.157 


3.03 


0.503 


-2.37 


5 


2.5 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.25 


1.36 


-8.08 


-3.06 


3.01 


0.596 


-2.41 


5 


2.5 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.23 


1.32 


-9.8 


-1.79 


3.06 


0.761 


-2.25 


6.5 


10 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


11 


6 


450 


3.94 


0.499 


2.24 


-8.27 


-3.85 


2.56 


0.762 


-0.746 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2.3 


35 


4 


2 


450 


1.86 


0.563 


0.548 


-5.88 


-0.194 


3 


0.775 


-1.54 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


6 


250 


3.64 


2.17 


0.784 


-5.92 


-0.357 


2.35 


0.869 


-0.185 


5 


10 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


0.972 


1.25 


-8.12 


-3.61 


3.2 


0.877 


-3.98 


5 


2.5 


0.4 


3 


0.65 


3.6 


1.7 


2.3 


35 


4 


6 


450 


3.76 


0.69 


0.496 


-5.98 


-0.18 


3.12 


0.934 


-0.987 


5.75 


2.5 


0.34 


1 


0.65 


3.6 


1.82 


2.54 


36 


4 


3 


450 


3.34 


0.539 


0.544 


-6.02 


-0.141 


3.09 


0.939 


-2.37 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


20 


7 


2 


250 


3.28 


0.994 


2.03 


-5.68 


-0.692 


1.54 


0.962 


0.137 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


11 


10 


450 


3.35 


0.922 


2.35 


-5.82 


-1.22 


1.8 


0.974 


0.0869 


8 


10 


0.45 


3 


0.65 


3.6 


1.8 


2.18 


35 


14 


3 


450 


2.68 


0.33 


1.2 


-6.13 


-0.136 


2.94 


1 


-2.59 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


35 


11 


10 


450 


3.49 


0.915 


0.753 


-5.66 


-0.363 


2.74 


1.04 


-0.505 


8 


2.5 


0.45 


1 


0.65 


3.6 


1.8 


2.36 


40 


14 


8.5 


450 


3.97 


0.521 


0.636 


-6.16 


-0.131 


3.14 


1.07 


-1.43 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


7 


6 


450 


1.92 


0.776 


1.52 


-5.71 


-0.456 


2.53 


1.08 


0.0529 


8 


2.5 


0.45 


1 


0.65 


3.6 


1.8 


2.36 


34 


14 


3 


450 


2.59 


0.495 


0.607 


-6.3 


-0.0546 


3.14 


1.1 


-1.75 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


6 


250 


1.71 


2.07 


1.72 


-6.76 


-3.36 


1.43 


1.14 


0.247 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


5 


450 


1.83 


0.863 


1.09 


-5.76 


-0.325 


2.26 


1.15 


-0.114 


6.5 


2.5 


0.4 


3 


0.65 


3.6 


1.7 


2.6 


35 


7 


6 


450 


2.02 


0.389 


0.386 


-5.76 


-0.238 


3.03 


1.16 


-2.52 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.6 


20 


7 


6 


250 


1.78 


0.965 


0.815 


-5.82 


-0.28 


2.74 


1.17 


-0.792 


5.75 


10 


0.34 


3 


0.65 


3.6 


1.82 


2.54 


36 


4 


3 


450 


3.34 


0.604 


0.346 


-6.21 


-0.0893 


3.26 


1.17 


-2.93 


5 


2.5 


0.4 


1 


0.65 


3.6 


2.15 


2.3 


35 


4 


10 


450 


3.49 


0.614 


0.505 


-5.96 


-0.168 


3.03 


1.19 


-2.28 


8 


2.5 


0.45 


1 


0.65 


3.6 


1.8 


2.18 


40 


14 


8 


450 


3.34 


0.609 


0.376 


-6.1 


-0.114 


3.16 


1.22 


-1.77 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


2 


250 


3.03 


1.81 


1.62 


-6.85 


-3.19 


2.39 


1.24 


0.0627 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


40 


14 


8 


450 


3.34 


1.33 


1.85 


-6.18 


-1.86 


2.36 


1.28 


0.194 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


6 


250 


3.8 


1.12 


2.37 


-6.26 


-2.61 


1.3 


1.28 


0.325 


8 


5 


0.4 


2 


1.6 


2.5 


1.8 


2.36 


40 


14 


8 


450 


3.24 


1.55 


1.14 


-6.45 


-0.204 


2.72 


1.28 


-0.423 


5 


2.5 


0.4 


1 


0.65 


3.6 


2.15 


2 


35 


4 


6 


450 


2.07 


0.648 


0.942 


-6.16 


-0.221 


3.08 


1.28 


-1.83 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


4 


6 


450 


3.97 


0.793 


1.93 


-9.05 


-3.33 


2.67 


1.29 


-0.0918 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


35 


7 


6 


450 


2.02 


0.635 


0.699 


-5.73 


-0.297 


2.84 


1.3 


-0.716 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


2 


250 


2.15 


2.26 


1.36 


-8.13 


-2.91 


2.25 


1.3 


0.354 


8.5 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


3 


450 


1.8 


0.875 


1.08 


-5.98 


-0.267 


2.75 


1.33 


-0.998 


5.2 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2.5 


450 


1.83 


0.804 


1.24 


-5.74 


-0.387 


2.7 


1.35 


-0.487 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


3.52 


1.51 


0.996 


-5.81 


-0.365 


2.51 


1.35 


0.113 


5.5 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


2.25 


0.795 


1.25 


-5.73 


-0.403 


2.74 


1.35 


-0.167 


5.2 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


4 


2 


450 


2.46 


0.972 


1.18 


-5.75 


-0.344 


2.27 


1.38 


-0.194 


5 


5 


0.36 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2 


450 


1.9 


0.86 


1.1 


-5.68 


-0.375 


2.5 


1.41 


-0.575 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.15 


35 


4 


2.5 


450 


1.63 


1.46 


1.83 


-6.27 


-2.25 


2.61 


1.42 


-0.0187 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


2.5 


450 


3.17 


1.05 


1.15 


-5.74 


-0.42 


2.77 


1.42 


0.195 


6.5 


2.5 


0.4 


1 


0.65 


3.6 


2.15 


2 


35 


11 


10 


450 


1.9 


0.636 


1.65 


-7.12 


-0.707 


2.96 


1.42 


-3.25 


5 


5 


0.36 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2.5 


450 


1.8 


0.974 


1.3 


-5.92 


-0.34 


2.76 


1.42 


-0.214 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


2.4 


450 


1.65 


0.951 


1.13 


-5.69 


-0.392 


2.6 


1.44 


-0.815 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


40 


14 


5.5 


450 


2.13 


0.627 


0.885 


-5.69 


-0.327 


2.59 


1.44 


-0.513 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


6 


250 


2.79 


2.3 


1.03 


-7.24 


-2.22 


2.65 


1.47 


0.57 


4.8 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


4 


2 


450 


1.87 


0.895 


0.975 


-5.69 


-0.344 


2.45 


1.48 


-0.488 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


2 


450 


1.69 


0.68 


1.32 


-5.69 


-0.448 


2.79 


1.49 


-1.58 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.1 


35 


4 


2 


450 


1.92 


1.55 


1.9 


-6.12 


-1.09 


2.19 


1.49 


0.0385 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


35 


14 


3 


450 


2.68 


1.02 


1.87 


-6.24 


-2.74 


2.73 


1.5 


0.373 


8.5 


2.5 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


7 


2 


450 


3.62 


0.973 


1.27 


-7.49 


-1.6 


3.08 


1.53 


-3.54 



TABLE XI: Same as for Table llXl but now for our gaugino-dominated DM candidate point. 
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Propagation Parameter Combination 



Burkert Profile 



D 


Bo 


a 


80 


e 90 


e si 


U 31 


n 92 




L 


dV 

dz 


4>@ 


i 

Xlar 


N 


M 


logio(A) 


7 


logio(BF) 


■2 

Xdm 


<5xdm 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2.3 


35 


4 


2 


450 


1.86 


0.599 


0.65 


-5.94 


-0.197 


2.04 


0.542 


-1.77 


5 


2.5 


0.4 


3 


0.65 


3.6 


1.7 


2 


35 


4 


6 


450 


3.97 


0.85 


1.63 


-6.84 


-0.106 


1.87 


0.587 


-0.993 


8 


10 


0.45 


1 


0.65 


3.6 


1.8 


2.18 


40 


14 


8 


450 


3.34 


0.417 


1.21 


-6.13 


-0.145 


1.87 


0.629 


-1.88 


5.75 


2.5 


0.34 


3 


0.65 


3.6 


1.82 


2.36 


38 


4 


3 


450 


2.92 


0.575 


0.743 


-5.87 


-0.236 


1.95 


0.691 


-1.32 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2.3 


35 


7 


10 


450 


3.78 


0.672 


0.72 


-6.26 


-0.0915 


2.07 


0.833 


-1.13 


5.75 


2.5 


0.34 


3 


0.65 


3.6 


1.82 


2.36 


36 


4 


2.5 


450 


3.17 


0.643 


0.661 


-5.96 


-0.181 


2 


0.871 


-1.41 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


4 


6 


450 


3.97 


0.893 


0.838 


-6.22 


-0.121 


1.93 


0.871 


-0.515 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


7 


10 


450 


3.78 


1.01 


1.37 


-5.7 


-0.431 


1.3 


0.877 


-0.186 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


2 


250 


2.62 


1.74 


0.88 


-6.08 


-0.251 


1.41 


0.879 


-0.321 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


6 


250 


3.64 


2.3 


0.979 


-6.51 


-0.123 


1.24 


0.91 


-0.144 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


7 


2 


250 


3.2 


1.94 


0.84 


-6.19 


-0.206 


1.31 


0.913 


-0.228 


5 


2.5 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.27 


1.34 


-7.82 


-1.22 


1.97 


0.915 


-2.09 


5.5 


10 


0.34 


3 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


2.25 


0.532 


0.665 


-5.95 


-0.212 


2.1 


0.942 


-1.47 


6.5 


10 


0.4 


1 


0.65 


3.6 


2.15 


2 


35 


7 


10 


450 


2.59 


0.52 


1.7 


-6.55 


-0.177 


1.85 


0.95 


-1.91 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


6 


250 


2.79 


1.35 


1.3 


-5.88 


-0.376 


0.869 


0.957 


0.0602 


8 


5 


0.4 


2 


1.6 


2.5 


1.8 


2.36 


40 


14 


8 


450 


3.24 


0.863 


1.17 


-5.69 


-0.385 


1.37 


0.971 


-0.107 


5.75 


2.5 


0.34 


3 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


3.52 


0.63 


0.764 


-5.94 


-0.183 


1.8 


0.979 


-0.63 


5 


2.5 


0.6 


3 


0.65 


3.6 


1.7 


2 


35 


11 


10 


250 


1.62 


0.534 


0.966 


-6.24 


-0.0951 


1.96 


0.982 


-3.11 


5 


2.5 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.26 


1.33 


-8.87 


-2.97 


1.97 


0.992 


-2.02 


5 


2.5 


0.4 


1 


0.65 


3.6 


2.15 


2 


35 


4 


6 


450 


2.07 


0.709 


0.886 


-6.28 


-0.107 


2.06 


0.997 


-2.11 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


42 


14 


8 


450 


2.79 


0.644 


0.926 


-5.63 


-0.361 


1.57 


1.11 


-0.646 


6.5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


2 


100 


3.81 


0.793 


0.564 


-5.81 


-0.271 


1.85 


1.11 


-1.2 


8.5 


2.5 


0.4 


1 


0.65 


3.6 


2.15 


2 


35 


11 


6 


450 


1.97 


0.542 


1.83 


-6.82 


-0.226 


1.83 


1.11 


-0.921 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


40 


14 


8 


450 


3.34 


1.35 


1.89 


-6.23 


-2.04 


1.22 


1.15 


0.0709 


8 


10 


0.45 


3 


0.65 


3.6 


1.8 


2.36 


34 


14 


3 


450 


2.59 


0.418 


0.495 


-6.05 


-0.118 


2.1 


1.16 


-2.63 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


6 


250 


1.71 


1.57 


2.01 


-6.8 


-3.18 


0.939 


1.2 


0.303 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


35 


11 


10 


450 


3.49 


1.44 


0.894 


-6.29 


-0.119 


1.69 


1.23 


-0.323 


5 


10 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.03 


1.11 


-7.73 


-1.96 


2.23 


1.24 


-3.61 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


20 


7 


2 


250 


3.28 


1.12 


1.92 


-5.7 


-0.707 


0.86 


1.25 


0.422 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


2 


250 


2.38 


1.32 


1.92 


-5.75 


-0.812 


0.812 


1.25 


0.259 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


3 


450 


2.93 


0.743 


0.974 


-5.76 


-0.323 


1.61 


1.25 


-0.836 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


11 


10 


450 


3.35 


2.86 


0.979 


-6.92 


-2.26 


0.754 


1.27 


0.387 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


2 


250 


3.03 


2.04 


1.5 


-7.19 


-2.56 


1.18 


1.29 


0.113 


5.5 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


2.25 


0.835 


1.19 


-5.67 


-0.41 


1.61 


1.29 


-0.268 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


3.5 


450 


1.68 


0.751 


1.33 


-5.77 


-0.344 


1.54 


1.31 


-0.519 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


2 


250 


2.15 


1.06 


1.81 


-6.12 


-0.443 


1.46 


1.31 


0.364 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


5.5 


450 


2.14 


1.46 


1.07 


-6.33 


-0.164 


1.58 


1.32 


0.232 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


34.5 


4 


2.5 


450 


1.67 


0.69 


1.34 


-5.63 


-0.438 


1.58 


1.33 


-0.65 


5.2 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


4 


2 


450 


2.46 


0.974 


1.14 


-5.73 


-0.366 


1.58 


1.33 


-0.244 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


3 


450 


1.64 


0.713 


1.16 


-5.65 


-0.365 


1.36 


1.34 


-0.231 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


4 


450 


1.9 


1.07 


0.879 


-5.8 


-0.319 


1.52 


1.34 


-0.561 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


4 


6 


450 


3.76 


1.4 


1.44 


-6.11 


-0.319 


1.5 


1.36 


-0.422 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


4 


2 


450 


1.98 


0.86 


0.708 


-6.41 


-0.0764 


2.05 


1.37 


-0.828 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


4 


2.5 


450 


1.96 


1.03 


1.05 


-5.73 


-0.34 


1.46 


1.39 


-0.378 


5.2 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2.5 


450 


1.83 


0.757 


1.31 


-5.75 


-0.36 


1.55 


1.4 


-0.442 


5.75 


10 


0.34 


1 


0.65 


3.6 


1.82 


2.18 


36 


4 


3 


450 


3.69 


0.721 


0.777 


-6.3 


-0.147 


2.17 


1.4 


-0.526 


5 


2.5 


0.4 


3 


0.65 


3.6 


2.15 


2 


35 


4 


10 


450 


3.8 


0.677 


1.66 


-7.89 


-1.54 


2.06 


1.42 


-1.1 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.15 


35 


4 


2.5 


450 


1.63 


1.4 


1.94 


-6.37 


-3.47 


1.52 


1.44 


-0.053 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.23 


35 


4 


2.5 


450 


1.63 


0.936 


1.41 


-5.81 


-0.405 


1.63 


1.45 


-0.28 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


4 


3 


450 


2.1 


0.875 


1.28 


-5.79 


-0.358 


1.58 


1.45 


-0.153 



TABLE XII: Same as for Table lVlTll but now for 



our mixed DM candidate point. 
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Propagation Parameter Combination 



Einasto Profile 



D 


Bo 


a 


«o 


e 90 


e 91 


n 91 


1lg 2 


"A 


L 


dV 

Hz 


<t>Q 


'2 

Xlar 


N 


M 


logio(^) 


7 


login (BF) 


'2 

Xdm 


SXdm 


5.75 


2.5 


0.34 


3 


0.65 


3.6 


1.82 


2.18 


36 


4 


3 


450 


3.69 


0.752 


1.3 


-6.45 


-0.0536 


1.38 


0.421 


-0.946 


8 


2.5 


0.4 


1 


0.65 


3.6 


1.8 


2.36 


40 


14 


8 


450 


3.24 


0.536 


0.664 


-6.06 


-0.138 


1.49 


0.747 


-1.41 


5 


2.5 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.26 


1.35 


-7.47 


-0.381 


1.39 


0.943 


-2.06 


5 


2.5 


0.2 


3 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.28 


1.22 


-7.73 


-0.517 


1.5 


0.944 


-2.06 


5.75 


10 


0.34 


1 


0.65 


3.6 


1.82 


2.18 


36 


4 


3 


450 


3.69 


0.691 


0.927 


-6.26 


-0.107 


1.39 


0.969 


-0.955 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


20 


11 


2 


250 


2.62 


1.75 


0.914 


-6.22 


-0.17 


0.871 


0.98 


-0.22 


5 


10 


0.2 


1 


0.65 


3.6 


2.6 


2 


20 


4 


6 


450 


3.67 


1.01 


1.2 


-7.83 


-0.659 


1.61 


0.996 


-3.02 


5.75 


2.5 


0.34 


1 


0.65 


3.6 


1.82 


2.54 


36 


4 


3 


450 


3.34 


0.507 


0.535 


-5.81 


-0.235 


1.39 


1.03 


-2.28 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.6 


35 


11 


10 


450 


3.49 


0.745 


0.792 


-5.6 


-0.367 


0.997 


1.03 


-0.518 


8 


10 


0.45 


3 


0.65 


3.6 


1.8 


2.54 


35 


14 


3 


450 


3.23 


0.349 


0.391 


-5.96 


-0.173 


1.62 


1.08 


-5.56 


6.5 


2.5 


0.4 


3 


0.65 


3.6 


2.15 


2 


35 


7 


10 


450 


2.59 


0.723 


0.968 


-6.38 


-0.106 


1.53 


1.13 


-2.66 


8 


10 


0.45 


3 


0.65 


3.6 


1.8 


2.36 


35 


14 


5 


450 


1.83 


0.45 


0.599 


-6.13 


-0.138 


1.65 


1.14 


-2.76 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


40 


14 


8 


450 


3.34 


1.02 


1.66 


-5.78 


-0.537 


0.696 


1.16 


0.0664 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


3.52 


0.782 


1.3 


-5.68 


-0.381 


0.865 


1.17 


-0.196 


5.5 


2.5 


0.34 


1 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


2.25 


0.621 


0.657 


-5.87 


-0.251 


1.41 


1.19 


-1.21 


5.75 


2.5 


0.35 


1 


0.65 


3.6 


1.82 


2.36 


36 


4 


3 


450 


3.77 


0.546 


0.67 


-5.76 


-0.293 


1.41 


1.21 


-1.14 


5 


10 


0.4 


1 


0.65 


3.6 


2.15 


2 


35 


4 


10 


450 


3.8 


0.703 


0.735 


-6.28 


-0.132 


1.62 


1.21 


-2.18 


8.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


11 


2 


450 


2.08 


1.48 


1.35 


-6.52 


-0.146 


0.827 


1.24 


-0.285 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


5 


450 


1.83 


0.844 


1.01 


-5.69 


-0.371 


0.991 


1.25 


-0.0153 


5.5 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


2.25 


1.07 


1.18 


-5.8 


-0.358 


1.05 


1.26 


-0.297 


6.5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


11 


10 


450 


3.35 


1.79 


1.73 


-6.1 


-1.67 


0.371 


1.29 


0.331 


5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


6 


250 


1.71 


1.63 


1.97 


-7.56 


-0.329 


0.29 


1.3 


0.405 


8 


10 


0.45 


3 


0.65 


3.6 


1.8 


2.36 


35 


10 


3 


450 


1.76 


0.401 


0.565 


-6.11 


-0.135 


1.63 


1.3 


-5.62 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.36 


35 


14 


4 


450 


1.9 


0.865 


1.09 


-5.9 


-0.28 


1.05 


1.31 


-0.516 


6.5 


2.5 


0.4 


1 


1.9 


2.95 


1.7 


2.6 


35 


11 


10 


450 


3.49 


0.538 


0.928 


-6.22 


-0.0806 


1.4 


1.32 


-2.94 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


11 


2 


250 


3.03 


1.95 


1.54 


-7.26 


-2.12 


0.747 


1.32 


0.149 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.15 


35 


4 


2.5 


450 


1.63 


0.937 


1.91 


-5.87 


-0.544 


0.831 


1.34 


-0.104 


8.5 


2.5 


0.4 


1 


0.65 


3.6 


2.15 


2 


35 


11 


6 


450 


1.97 


0.555 


1.76 


-7.51 


-0.198 


1.4 


1.37 


-0.666 


6.5 


5 


0.4 


2 


1.6 


2.5 


2.15 


2.3 


20 


7 


2 


250 


2.15 


1.73 


1.66 


-7.84 


-0.648 


0.854 


1.38 


0.429 


5 


10 


0.4 


3 


0.65 


3.6 


1.7 


2 


35 


7 


10 


450 


3.21 


0.597 


2.02 


-8.11 


-1.95 


0.98 


1.39 


-0.174 


8 


10 


0.45 


1 


0.65 


3.6 


1.8 


2.36 


35 


14 


4 


450 


1.9 


0.421 


0.521 


-6.04 


-0.128 


1.47 


1.46 


-2.32 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


30 


4 


2 


250 


3.36 


1.34 


1.72 


-6.42 


-0.22 


0.766 


1.47 


-0.146 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.54 


36 


4 


3 


450 


3.34 


0.745 


0.784 


-5.67 


-0.353 


1.27 


1.47 


-0.688 


5 


5 


0.38 


2 


1.6 


2.5 


1.81 


2.25 


35 


4 


2.5 


450 


1.65 


0.841 


1.24 


-5.71 


-0.385 


1 


1.47 


-0.607 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


20 


11 


6 


250 


2.79 


2.07 


1.23 


-6.84 


-3.26 


1.03 


1.48 


0.582 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


35 


7 


10 


450 


3.78 


1.75 


1.34 


-6.7 


-0.244 


1.11 


1.5 


0.441 


5 


5 


0.4 


2 


1.6 


2.5 


1.7 


2.3 


40 


4 


2 


450 


2.01 


0.58 


0.968 


-5.65 


-0.356 


1 


1.51 


-0.7 


4.8 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.3 


35 


4 


2.5 


450 


2.19 


0.841 


1.03 


-5.74 


-0.34 


1.14 


1.52 


-0.569 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


2 


450 


1.69 


0.78 


1.37 


-5.76 


-0.419 


1.15 


1.52 


-0.208 


5 


2.5 


0.4 


1 


0.65 


3.6 


1.7 


2 


35 


4 


6 


450 


3.97 


0.949 


0.769 


-6.43 


-0.0506 


1.37 


1.52 


0.137 


5 


10 


0.4 


1 


0.65 


3.6 


2.15 


2.3 


35 


4 


6 


450 


1.88 


0.5 


0.444 


-6.12 


-0.1 


1.58 


1.52 


-3.38 


8 


5 


0.45 


2 


1.6 


2.5 


1.8 


2.18 


35 


14 


3 


450 


2.68 


1.16 


1.81 


-6.38 


-3.26 


1.08 


1.53 


0.409 


5 


5 


0.38 


2 


1.6 


2.5 


1.8 


2.25 


35 


4 


3 


450 


1.65 


1.01 


1.28 


-5.79 


-0.407 


1.11 


1.54 


-0.938 


5.75 


2.5 


0.34 


3 


0.65 


3.6 


1.82 


2.36 


36 


4 


2.5 


450 


3.17 


0.499 


1.26 


-6.3 


-0.0601 


1.03 


1.54 


-0.744 


5 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


3 


450 


2.6 


0.932 


0.928 


-5.67 


-0.406 


1.26 


1.55 


-0.258 


5.75 


5 


0.34 


2 


1.6 


2.5 


1.82 


2.36 


36 


4 


2.5 


450 


3.17 


1.18 


1.16 


-5.78 


-0.436 


1.18 


1.55 


0.179 


8 


10 


0.45 


3 


0.65 


3.6 


1.8 


2.36 


40 


14 


8.5 


450 


3.97 


0.495 


0.295 


-6.01 


-0.134 


1.63 


1.56 


-1.63 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2 


35 


11 


10 


250 


1.62 


0.965 


2.11 


-6.05 


-0.942 


0.881 


1.58 


0.0818 


5 


5 


0.6 


2 


1.6 


2.5 


1.7 


2.3 


35 


11 


10 


250 


1.53 


0.477 


0.935 


-5.66 


-0.377 


1.32 


1.59 


-0.686 


8.5 


10 


0.4 


3 


0.65 


3.6 


2.15 


2 


35 


11 


6 


450 


1.97 


0.484 


1.68 


-7.26 


-1.91 


1.54 


1.6 


-1.12 



TABLE XIII: Same as for Table llXl but now for our mixed DM candidate point. 



